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ABSTRA CT 

A thyristor valve for HVDC syfeems must contain a large 
number of devices connected in series and in parallel to 
accomriodate large voltage . and current requirements.' As the 
turn on and turn-off char act eri tics of the individual thyristors 
vary, the thyristors are subjected to overvoltages and over- 
currents during the turn-on. and commutation of a valve. The 
overvoltages are controlled by the grading and damping circuits. 
Computer programs are developed to calculate turn-on and turn- 
off overvoltages in a J-IVDC valve having thyristors connected 
in series. 

Power loss in the damping circuit is a significant 
component of the overall convertor losses and hence has a 
bearing on its design. A method is developed for the exact 
calculation of losses in the damping circuit. The 'effect of 
bridge operation and damping circuit parameters on losses is 
analysed. 

A comptiter program is developed to calculate the 
sensitivities of the valve voltage v/ith respect, to damping 
circuit resistance and caacifance in Time Domain using Adjoint 
netxTOrk approach'.' A metho'"" of parameter optimization i.s used 
to select the values of damping circuit parameters which reduce 
the recovery voltage without undue increase in losses. 



CHAPTER 1 


INTRODUCTION 


Status of HVTX: Transmission 

Despite the general acceptance of A.C, transmission^ 
D.C’i' Transmission is preferred for certain applications 
because of its advantages. In particular it can be consi- 
dered for the following applications. 

(1) Submarine cables of length more than 20 miles, 

(2) Interconnecting A.C. systems of different 
frequencies , 

( 3 ) For transmitting large amounts of power over 
long distance by overhead lines. 

( 4 ) For transmission by cables in congested 
urban areasV 

Ratings of HVDC link have grown continuously and at 
present a power rating of more than 2000 MW with transmission 
voltage upto + 600 KV are feasible. The use of better 
converter transformers^ harmonic filters, sophisticated 
methods of converter control and protection, 12 pulse 
converters has improved the performance of HVDC to a great 
extent. The modem trends towards the use of multiterminal 
systems and the developments in the design of HVDC circuit 
breakers will certainly increase the application of HVDC 


in futiarel' 



2 


1.2 HVDC Thyristor Valve [7, 13-16] 

The major component of the converter station is the 
valve which allows current flow only in one direction and 
the instant of starting the conduction through a valve can 
be controlled with the help of grid (firing angle) control. 

The very high current density of the cathode spot in 
mercury arc valves / together with high breakdown level of 
me%3ury vapour, established the superiority of mercury arc 
valves during early thirties of this century. Mercury 
arc valves are now displaced by thyristor valves which 
have proved to ' be more .reliable. In particular, with 
thyristor valves 

1) there is no problem of arc backs (reverse conduction) 

2) the converter is of modular construction and the 
design and maintenance are easy 

3 ) the operation is more reliable and capable of fast 
controls 

A thyristor valve consists of a number of individual 
thyristors connected in series and parallel to accommodate 
the large voltage and current requirements because of the 
llndted voltage and current carrying capability of 
individual thyristors. The thyristor is a three terminal, 
semiconductor device made of alternate layers of p and 
n type silicon'.' The diagram of the device and electric 
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symbol Is shown in Fig.1.1",’ The device can carry current 
only in one direction from anode to cathode and the instant 
of initiation of conduction can be controlled by gate".' 

The characteristic of the device is shown in Fig, 1,2 
with the three states of the device (l) Reverse bias and 
blocking (2) forv/ard bias and blocking and (3) forward bias 
and conducting,' 

A thyristor is turned on by applying an electrical 
signal to gate when the device is forward biased. Voltage 
and current waveforms during S.C.R, turn-on time are shown 
in Fig,l','3. t'^hen firing current is delivered to the gate of the 
device. It does not switch immediately from the blocking to the 
full conduction state. Indeed for a short period of time, 
the device continues to block the anode voltage applied to 
it in almost the same way as it would had the firing pulse 
net yet been delivered. Thereafter, the forward impedance 
starts to decrease, but not until a further time period has 

elapsed does the device became fully turned on. The total 

1 . 

turn-on time of the device is subdivided into two distinct 
periods, called the delay time and the rise timeT Both the 
delay time and the rise time are related to the rise time and 
amplitude of the gate firing current. The total turn-on time 
is of the order of 2-4 microseconds',,' 

Turn-off of a device implies that the forward conduction 
has stopped and the application of a positive voltage to the anode 




LHC Tk;\C/> (_ O7O160L cr Ti*7kISTo< 


4 
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T-itVJClSTof^, 





5 


will not cause any current to flow if there is no gate 
signal,' Voltage and current waveforms during SCR turn-off 
time are shown in Fig'., 1.4, The total turn-off time is 
subdivided into two regions , called the reverse recovery 
time, and the gate recovery time. When a reverse voltage 
is applied across the thyristor its current reduces to zero 
value. Because of the storage charge present in the device, 
it continues to conduct in the reverse direction until its 
storage charge is neutralised. The reverse recovery time 
is a function of forward current and the rate of decay of 
forward current . 

During the gate recovery time, reverse voltage 
must be maintained .across the device'. Not until the end 
of the gate recovery time is the device capable of blocking 
reapplied forward anode voltage'^ Th:|.s- time is dependent 
upon Junction temperature of the devj.ee and the rate of 
reapplication of forward voltage'. . 

The current rating of a thyristor depends upon the 
dimensions of the device as well as cooling arrangements. 

The voltage rating of a thyristor indicates the maximum 
reverse voltage it can block v^ithout excessive reverse 
current'. At present, the thyristors used in HVDC valves 
J^ave a rpaximum voltage rating of 3,5 KV and a current 
rating qf leOO Amps. 
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The proper sharing of voltage by the Individual 
devices in the series connected thyristors in steady state 
is achieved by connecting a resistance or a zener diode 
across each device. Transient voltage sharing is ensured 
by a grading circuit which consists of series connected 
resistance and, capacitance connecting across each device. 

Cooling of a HVDC valve is required because the 
thyristors are thermally limited devices, A variety of 
fluids such as air^ oil, water and SF^ can be used for this 
purpose. Air cooling is simpler and more reliable. The 
insulation of the valve can also be through air, oil or 
SFg gasV 

Overvoltages in a HVPC Valve 

The electrical design of HVDC valves depends upon 
the overvoltages generated in the system and those generated 
within the valve, The valve cS subjected to internal over- 
voltages and voltage stresses caused by the external system 
conditions’. Switching surges from the A.C’, side and voltage 
surges from D,C, side may cause overvoltages across a valve. 
The overvoltages can also occur during abnormal conditions 
(maloperation of the valve) such as arc back, arc through, 
quenching and misfire".' 

The turn-on and turn-off overvoltages constitute 
internal overvoltages in a valve",' As the characteristics 
of thyristors vary due to manufacturing tolerances, the 
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individual thyristors turn-on and turn-off at different 
times because of v;hich the thyristors and the valve are 
subjected to overvoltages, .hnd also, a recovery overvoltage 
is developed across the valve after comi.utation. Grading 
and damping circuits are used to limit those overvoltages. 

1.3 Obj ectiv e of ^ The s is 

The objective of this thesis is to study in detail 
the turn-on and turn-off overvoltages in a HVDC Thyristor 
Valve and design of damping circuit. The turn-on and turn- 
off processes in a valve have to be simulated accurately 
to Icnov' the overvoltages developed across it. Digital 
simulation/ based on nepwork equations derived from topo- 
logical considfirations / is used to calculate these over- 
voltages. Effect of the valve and the system paramctv'.rs is 
also investigated.' 

The considerations in the damping circuit design are 
overvoltage and power loss. An accurate evaluation of the 
ov'rvoltaye and power loss and the variation with damping 

circuit parameters is essential in the optimal de- ign of 
the damping circuit. A rnethematical technique to optimize 
the paramot'')rs is also investigated, 

RGview of th e ^ L itera ture 

A detailed method of calculating turn-on overvoltages 
in a HVDC valve having thyristors in series is presented in 
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refV[l] by G. Karady and T. Gilsig/ in v;hich a cumulative 
delay time distribution curve is used to select the number 
of thyristors that are turned on at any instant’. The thyristor 
is modelled by a voltage source in series with a switch v^hich 
closes; at .the instant of delay tim.e of the thyristor’.’ During 
the rise time the voltage across a thyristor will decrease 
exponentially’.’ The effect of various system and valve para- 
meters on turn-on overvoltages is analysed taking a represen- 
tative thyristor valve. An IBM supplied ECAP pro^^etm is used 
for circuit analysis. An extension of this work is carried 
out in ref .[2],, in vj-hich parallel thyristor strings are also 
considered, which, are necessary to moot the current require- 
ments. Effects on overvoltages, of parallel strings of 
various turn-on delay time distributions and the effect of 
series string inductances on current sharing and overvoltages 
are considered’. The system is described by a sot of state 
equations and solved numerically. 

Ref, [3] gives an exact method of calculating turn-off 
overvoltages which consist of both the overvoltages caused due 
to non-simultaneous turn-off of individual thyristors and the 
recovery overvoltages occurring during and after the turn-off 
process’," A statistical distribution of storage charge is 
assumed for the thyristors’." The thyristors are divided into 
a number of groups and each group of thyristors is represented 
by a step change in impedance from a small value to a very high 
value'.’ The results are evaluated using ECAP program. 
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In ref. [ 7 ], a simplified analysis is carried 
calculate the valve voltage waveform. In ref .[^4]/ a detailed 
equivalent circuit is presented to calculate the recovery 
voltage, neglecting the turn-off process. 

Turn-on and turn-off overvoltages can be controlled 
by providing proper grading and damping circuits. Because 
of the presence of voltage across a vulve during its non- 
conducting state the losses take place in the damping circuit. 
An approximate formulae are suggested by Ainsworth [17] to 
calculate the valve damping loss. An improvement to this 
method is made in ref, [4] , by taking a second order model 
for generating voltage oscillations, which are superimposed 
on the ideal valve volt age V A hybrid computer is used to 
evaluate the power loss',' 

The design of the damping circuit requires the 
determination of the capacitance and resistance which limit 
the recovery voltage to a predetermined level without undue 
increase in losses. The various considerations to bo taken 
in designing a damping circuit arc given in ref ,[4] .' 

V 

1,5 Summary of the Thesis 

■ A chaptorwise summary of the thesis is given below. 

Chapter 2 deals with the calculation of turn-on 
overvoltages and the effect of system and valve parameters 
on it7 A general computer program has been developed to 


1 
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simulate the turn-'On process of thyristors In a valve. The 
approach adopted for simulation of a valve during turn-on 
period- is along the linos given in refs'. [1,2]. The 
equivalent circuit is simulated by developing the state 
equations on the basis of tfcpological considerations. 

A method of calculation of turn-off overvoltages 
using a digital simulation is presented in Chapter 3 . 

During the turn-off period the valve is represented by a 
time varying resistance in series wi-th a time varying 
capacitance based on the approach given# in ref . [s] . 

The effect of system and valve parameters on turn-off 
overvoltages is studied'. The recovery voltage neglecting 
turn-off process is also calculated'.' 

Chapter 4 deals v/ith the calculation of sensitivities 
of the valve voltage with respect to damping circuit resis- 
tance and capacitance using adjoint network approach. , The 
effect of damping circuit parameters on the sensitivities 
has also been studied •' 

An exact method of calculating the damping circuit 
loss is presented in Chapter 5# using the equivalent circuit 
representation of the converter system during different 
intervals of conduction’.' The results are compared with the 
results obtained using an approximate second order system 
model for valve voltage oscillations and with the results 
obtained by the approximate formulas given by Ainsworth'.' 
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A method of evaluation of optimal damping circuit 
parameters is presented in Chapter 6 , by minimizing an 
objective function which takes into account both recovery 
voltage transient and damping circuit loss'.' 

Conclusions of this thesis are given in Chapter 7. 
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CHAPTER 2 

CALCULATION OF TURN-ON OVERVOLTAGES 

2'ri Introduction 

A thyristor valve used for HVDC transmission systems 
must contain a large number of thyristors connected in series., 
because of the limited voltage capabilities of a single 
thyristor. Because of the unequal turn-on characteristics 
of the individual thyristors/ the thyristors turn-on- non- 
simultaneously due to which the thyristors are subjected to 
ovebvoltages during this transitional period. A computer 
program is developed to calculate the turn-on overvoltages in 
a HVDC valve having number of thyristore in series. The effect 
of critical system and valve parameters on turn-on over- 
voltages is studied by considering a numerical example. This 
study is useful for valve design and for optimization of system 
and valve parameters'^ 

2.2 Equivalent Circuit for the Turn-On Process of a 
Thyristor Valve . 

2.2.1 Equivalent Circuit of the Supply Networks 

The thyristor valve converter system [l]/ is shown in 
Figi'2.1. Each thyristor valve will operate as one leg of the 
three phase bridge circuit. Each valve is shunted by the 
damping circuit/ which consists of a resistance (R^) and a 
capacitance (C )/ and has in series a parallel combination 
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of R -L / which is known as anode damper. The distributed 
Q a 

transformer and switchyard capacitances are treated as 

lumped elements C. and C , respectively. 

t y 

During the normal operation of the three phase bridge, 
tv;o and three valves conduct alternately. The curre^- 
conduction through valves 3 and 4, is shown by the solid line 
in Fig. 2.1. At -the end of the conduction period of valves 3 
and 4, valve 5 starts to conduct resulting in three valves 
(3,4 and 5) conduction during the commutation period'. The 
current transfer from valve 3 to valve 5 during comrautation 
is shown by the dotted line in Fig. 2.1. During this 
commutation period, the voltage across valve 5 reduces to 
zero, and the current through valve 5 increases to I^. The 
transient distribution of the collapsing forward voltage 
across the individual thyristors of valve 5 is studied in 
this chapter".' 

In practice, the turn-on period is a few microseconds 
which is small compared to the time constant of damping 
circuit'. So there is no effective connection between phase 
a and phases b and c. Accordingly, valve 5 turn-on process 
is represented by the dotted path of Fig. 2.1 with the 
associated stray capacitances. A further simplification is 
done, through a star-delta transformation of capacitances, 
and neglecting circuits across valves 3 and 5 as 

their effect can be neglected compared to circuits and 

valve stray capacitance (Cy)'. The stray capacitance across 
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the valve 5 is neglected* The simplified equivalent circuit 
of the supply network is shown in Fig *2.2. The two voltage 
sources v-^ and v^ are combined to an equivalent voltage 
source V / which is expressed as V = v, -v where v and v 
are b and c phase voltages ^ using Thevenin's theorem* Since 
the peak voltage across the valve is of interest in calculating 
turn~on overvoltages ^ is assumed to be at the maximum value 
corresponding to a = 90^, in this analysis. The delta equi- 
valent capacitances terms of the star 

connected capacitances and are given as 


"be 


= V<2C^ 




= Sc ' 


( 2 . 1 ) 


The capacitance C-^ which comes across valve 3 is neglected. 


2.2.2 Equivalent Circuit of a Thyristor Valves 


A thyristor valve having N number of series connected 
thyristors /With a voltage grading circuit consisting of 
and O’ across each thyristor is sho^vn in Fig. 2. 3. Prior to 
the turn-on of each thyristor/ a positive voltage is developed 
across each thyristor. Accordin gly /during turn-on process / 
each thyristor can be represented as a voltage source/ which 
starts decreasing as the gate signal is applied. The voltage 
across a thyristor does not start to collapse until a finite 
time after the gate pulse is applied/ because of the f/lnite 
speed of charge carriers/ and takes some further finite time 
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to collapse to the conduction value. Accordingly, turn-on time 
(t^j^) is the sum of delay time (t^) and rise time (t^) which 
is shown in Fig. 2. 4. 


The parameters t^ and t^ for a given type of thyristor, 
under specific test conditions and for a given gate current 
waveform can be obtained from the specifications given by the 
thyristor manuf ■'cturers. In this analysis, the thyristor is 
represented by a switch closing at the instant t^, in series vjith 
a voltage source given by 


where 



( 2 . 2 ) 


(2.3) 


For a given type of thyristor, there will be a range 
of valtti^s of delay time because of manufacturing tolerances. 
Therefore, a statistical distribution is considered for the 
delay times. A cumulative delay time distribution curve, 
based on the measurements of a sample taken from a particular 
thyristor type, is shown in Fig. 2.5 [l]. 


Using the curve shown in Fig. 2. 5, we can obtain the 
number of thyristdrs in the on condition as a function of 
time after the gate pulse ig applied. The curve is approxi- 
mated as a finite number of steps, where each step represents 
a group of thyristors and the time at which that group turns-on. 


i 
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For a given thyristor type, the rise times between 
individual thyristors also vary, which cause second order 
effect on turn*-on overvoltages. A representative value of 
t^ can be chosen to calculate the val'uo of 6 v/hich is used 
in Eq. (2,3). For a number of high povjer thyristors, 6 
ranges from 2.0 - 3.5 x 10^ sec“^ [l]V 


2,2,3 Combined Equivalent Circuits 

The complete ecfuivalent circuit for the supply system 
and thyristor -valve is shown in Fig. 2.6,' 


2.2.4 Simplification of the Equivalent Circuits 

The equivalent circuit shown in Fig. 2, 6 is simplified 
for computer calculation and is shown in Fig. 2. 7. During the 
turn-on process the valve is represented by a time varying 
voltage source in series with a time varying and C^l 

The total number of thyristors (N) are divided into 
n groups'," The i’th group with thyristors is represented 
by a svj-itch which closes at t ^^. , in series with a voltage 
source V^(t) = N^',v^(t), Using the switching function of 
Fig. 2. 5, the time t _ . for each group of thyristors is obtained. 


At t < i.e. prior to the turn-on of first group 

of thyristors, the voltage across each thyristor is v^= 

The equivalent circuit of the valve during the period 
is shown in Fig, 2 .8(a) Group 1 having turned-on, it is 
represented by a voltage source 


d2 



22 





L-Ci 

— 

(-S 




w:i 

.^-?r?ofv 

^--iO 


K'£^l 


Lh 


(Cbc ^x; 


'*■ 

4-- 


n V 



.,1 ■'^1 pUP ‘ir ELQO# !\)T CiRcuir 

calculat) ojNj 


,7 


FOIn 



23 




M, V , 

1 ol 


~ 6 ) 


(2-4) 


where ~ ^o* remainder of the valve is represented by 

(N-Nj^)Rg in series with C g/(N-Nj^) with an initial voltage 
of Vgj^Ct^^) where 




(2.5) 


At t = t^ 2 '' 2 turns-on and the equivalent circuit 

of the valve is shown in Fig- 2-8(b). Group 2 is represented 
by the voltage source 






(2.6) 


where v ^2 voltage across a non-conducting thyristor 

which is given by 


''02 = [^'*^32’ ''g <“-h> + '"gl<'^d2’h‘"-h> 


The remainder of the valve is represented by (N-N^+N 2 )Rg 


in series with a capacitor Cg/(N— N^+N 2 ) with an initial 
voltage of Vg2^^d2^ where 


\2^’^d2^ 


V 02 (N - N^+N2) 


( 2 . 8 ) 


The above process continues through the switching of successive 
groups at t^ / t^ 4 / ^d(n-l)* during the period 

t - ■ . < t < t , the equivalent circuit of the valve is 

d(n— 1) dn 

shov/n in Fig'i' 2.8(c)'. The voltage applied across last 
group of thyristors at t = Is obtained from this circuity 
As a measure to the overvoltage across each thyristor an 
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overvoltage factor F(t) is defined as 


F(t) = v^/v^ (2.9) 

where v^ is the voltage across each thyristor in ith group 
given as 


= [lift) Rg(N i V^^(t)]/(N- N^H-Nj+...+N._^) 

( 2 . 10 ) 

Vq is the initial voltage across each thyristor prior to thc 
turn-on of first group and the time t is, in the period 


t ,*r*r: < t < t , . . 
di-1 di 


213 Calculation of Turn-On O/ervoltages 

2.3'.i System State Equations; 

The state equations are derived for using the 

algorithm given in Appendix A, The final state variable 
equations are given by Eq. (A .13 ) of Appendix A. 

2. '3. 2 Computer Program; 

The state equations are solved numerically by using 
Runge-Kutta fourth order method. The flow chart of the 
comput'^r program is shown in Fig. 2V9V 

2 ."4 An Example , 

2.4.1 Description of the Problem; 

A valve which contains 200 thyristops and which 
operates at a peak repetitive voltage of 200 KV is conside.r .• .r 
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IT 


vSTART 




READ input data 


set first group 
to conduction 





Fig. 2. '9? Flow chart of the computer program 
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to show the effect of certain critical parameters on tum-on 
overvoltages with the following supply system and valve 
parameters'. This example is taJcen from ref,[l]. 

Supply system: 

V = 200 KV 

o 

2L^ = 0.15 H' 

= 2.0 nF . = 2.0 to 20.0 nP 

= 1 mH 

cL 

R =4 M-ohm. 

a. 


Valve ; 


R = 1 to 100 ohms 
<3 

C = 0*1 ■to 10 micro F. 

6 = 2.0 X 10^ sec“^ . 


Number of thyristor groups considered 
in a valve ' = 14',' 

The number of thyristors in each group and their switching 
times are shown in Table 2.1V 

The step length used in Runge-Kutta siibroutine is 
0.01 microseconds, and the initial state variable vector 
is as follows'. 

V 2 (o) = 200,000 volts 
- 200,000 volts 
V (o) = 200,000 volts 

g 

iyCo) = 0 

lg(0,) =0 , 
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TABLE 2.1s Number of Thyristors in each Group and Their 
Switching Times'," 


Group Number of thyristors Switching time 

in each group for each group 

(mi croseconds ) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


3 

0.2 

4 

0.4 

7 

0.6 

16 

0.8 

30 

1.0 

32 

1.2 

31 

1.4 

27 

1.6 

16 

1.8 

14 

2.0 

7 

2.2 

6 

2.4 

4 

2.6 

3 

2.75 


14 
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2 .4 V2 Results: 

The computed valve voltage and valve current waveforms 
during turn-on' process are shown in Fig. 2.10. , The valve 
voltage is reduced from 200 KV to zero, and the valve current 
has increased from zero to 291 Amps, during the turn-on 
process of the valve',' 

The effect of valve param-'^ters and on turn-on 
overvoltages : 

The variation of overvoltage factor with time for 

different R is shown in Fig. 2,11, For a particular R and 

g y 

C the overvoltage factor increases v/ith time, as the voltage 
across a late firing thyristor is the sum of the resistive 
and capacitive components . Fig. 2'ill shows that reducing the 
value of R can reduce the overvoltage factor significantly. 

g 

For a value of R =1 ohm, the overvoltage factor is much 

Q 

less compared to the higher values of R . But the value of 
Rg as low as 1 olim is impractical, as this would result in 
excessi've values of di/^i- during turn— on, and insufficient 
damping during transients [l] • Hence a practical minimum valuA 
is taken as of the order of 5 - 10 ohms. 

The variation of overvoltage factor with -time for 
aifforent is Shown In Pig. 2.12. The effect of Increasing 
the value of is not significant at lower times, since the 
capacitance voltage is proportional to the time integral of 







voltage factor 
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the current. For a value of R = 10 ohms, increasing C 

g g 

beyond about 1,0 micro farad is not useful/ since the 

dominant component of voltage is due to resistance. Increase 

in C beyond a certain value is not advisable/ since the 

y 

grading circuit losses increase with increase in [l]«' 

The effect of svjitchyard capacitance: 

The variation of overvoltage factor with time for two 

different values of C is shown in Fig, 2.13. Decreasing 

yc 

C from 20 nF to 2 nF reduces overvoltage factor and also 
yc 

the reduction in overvoltage factor increases with time. 

The effect of niaraber of thyristors (N) : 

The yariation of -overvoltage for different N is shown 
in Fig. 2, 14V The overvoltage factor is reducing with decreasing 
Ni This is an important observation which .can be followed in 
building a valve with the highest voltage thyristors available . 

2,4,3 Discussion: 

In ref,'[l]/ for the transient solution, an IBM-modified 
FCAP (Electronic Circuit Analysis Program) program' is used as 
a, .'.subroutineV A data preparation program written in FORTRAN 
and a monitor program to link this with ECAP program are also 
usedV A complete solution for a valve containing 200 thyristors 
divided into 16 groups took 276 seconds c.p.u". time on an 
IBM 360-50 computer [l];' 


4 



Over voltage factor 




0 .15 H 
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In this analysis/ the system is represented by a set of 
state variable equations and are solved numerically', A complete 
solution for a valve containing 200 thyristors divided into 
14 groups took 0.'66 seconds c.p.u. time on DEC~1090 computer 

The results obtained in this analysis agree with the 
results of ref'.fl] qualitatively/ but slightly differ 
quantitatively';' This difference is expected to be due to 
difference in the value' of 6 and number of groups considered 
in the two analyses";' 

2.5 Conclusions 

A general computer program to calculate the turn-on 

overvoltages in a HVDC thyristor valve has been developed, 

which simulates the tum-on process with the scatter in 

delay times of the series connected thyristors of the valve'. 

The effect of varying the valve parameters R and C for 

y g 

a representative valve has been demonstrated. Minimizing 
the resistance of the voltage grading circuit decreases 
thyristor overvoltages. Increasing the value of capacitance 
has the same effect/ but the incremental improvement is 
small for practical values of resistance';' 

The valve current during turn-on period is supplied 
mostly from stray capacitances, the transformer current 
being small. For t'ue times of the; order of 1-3 microsec end e, 
the switchyard capacitance will act as a principal current 
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source'.' The current supplied by the switchyard and trans 
former capacitances is limited by the parallel connected 
anode resistance and inductance, the latter being more 
effective^ 

Thyristors of the highest available voltage rating 
should be used, to reduce turn-on overvoltages. 
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CHAPTER 3 

CALCU L/vTION OF TU RN- OFF OVSRVOLT7\GES 

3.1 I n t_r qd u c 

A thyristor valve used for H\/DC systems must contain a 
large number of thyristors connected in series, because of the 
limited voltage capability of a single thyristor. Because of 
the difference in turn-off characteristics of the individual 
thyristors, the thyristors turn-off non-simultaneously, 
because of which they are subjected to overvoltages during 
this transitional period'. Also, the valve is subjected to 
recovery overvoltages, which occur during and after the turn-off 
period. A computer program is developed to calculate the turn- 
off overvoltages which consist of both recoveiry overvoltages 
and overvoltages caused due to non— simultaneous turn-off of 
thyristors in a HVDC valve. The effects of critical system and 
valve parameters on turn-off overvoltages . are studied by taking 
a representative thyristor valve. This study is useful for vaJv® 
design and for optimization of system and valve parameters, A 
method of calculating the recovery overvoltage neglecting turn- 
off process has also been presented. 

3.2 Equiv alen t Circu it for t he Tu rn-Off Process of A Vg-lvg. ^ 

3.2.1 Equivalent Circuit of the Supply Systems 

The thyristor valve convertor system is shown in Fig,3.1 
"si , Each valve is shunted by the damping circuit 
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has In series a parallel combination of R -L / wluch Is known 

a a 

as anode damper. The distributed transformer and sv'.dtchjard 
capacitances are represented by the lumped parameters and 
Cy, respectively. 

During the normal operation of the three phase bridge 
circuity two and three valves conduct alternately. The current 
conduction through valves 3 and 4 -is shown by the solid line 
in Fig. 3.1. As the valve 5 turns on, the current through valve 3 
decreases to zoro, while the current through valve 5 increases 
to as shown in Fig, 3. 2 (a). During this commutation period, 
the voltage across valve 5 decreases to zero in a few micro- 
s .Gonds, while a rev-rse voltage is developed across the valve 3 
as shown in Fig,3.2(b). The turn-off of valve 3 is considered , as 
a representative case to calculate the turn-off overvoltages. 

Even after the current through valve 3 goes to zero, 
valve 3 continues to conduct in the reverse direction, because 
of the presence of storage charge represented by unneutralised 
carriers in the thyristors of valve 3 . The reverse current of 
peak value which also flows as a forward current through 
valve 5 superimposed on is shown by the dotted line in 
Fig. 3.1'. The time during which reverse current flows is some 
tens of microseconds and. is enlarged for clarity in Fig. 3. 2. 

, * 

As the values of storage charge for individual thyristors 
3iffer, the thyristors of valve 3 turn-off sequentially, thus 
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Increasing the impedance of the outgoing valve at a finite 
rate. The voltage across valve 3 continues to increase till 
all thyristors are turned-off. The transient voltage across 
valve 3 v;ill oscillate initially as shown in Fig. 3, 2(b), and 
will finally approach the steady state line voltage 


A valve consists of a number of series connected 
thyristors with a voltage grading network R across each 

y y 

thyristor. The eguivalent circuit of valve 3 turning-off 
is shown in F'''g.3.3. It is developed from Fig .3,1^ by repre- 
senting each of the non-conducting valves (6, 2 , and 1) by its 
damping circuit parallel with its grading network 

R “C , and the conducting valves, 5 and 4, by short circuits. 

The anode damper R -L is neglected for all the valves / except 

a a 

for valve 3. The capacitances and of Fig. 3. 3 

are given by the follov;ing expressions obtained from star-delta 


cai>?.citance transformations 

Sec 


'ab 


- 

3Ct +nCy 


(3.1) 


(a + c )' 

P z... 

^ac 3C^ +' 2C 


(3 .2) 

The maximum turn-off overvoltages occur when the yol+'ag>. 


across the valve is maximum. This -corresponds to' the instant 
cc + = 90°/ and when v_^ = = '~E^/\/'2/ and v^ = +Ej^/V2o 

In Fig.3.3/ the voltage is assumed to be constant equal to 
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E„/'/2, as the variation in phase voltages is small during the 
L) 

turn-off process. The initial rate of change of current is 


given by 


dt L^" 


(3.3) 


3.2.2 Equivalent Circuit of a Thyristor Valve 


Even after the current through a thyristor falls to 
zero, it will continue to conduct, because of unneutralised 
carriers present in it. The negative current-time integral 
of Fig'.'3'.2 (a) is shown in detail in Fig. 3. 4. It is assumed 
that the thyristor turns-off at the instant when the reverse 
current reaches its peak, and a reverse voltage is suddenly 
applied across it. Accordingly, the storage charge is taken 
as the shaded area of Fig.3.4; and each thyristor is represen- 
ted as an ideal svyitch v/hich is opened at the instant when the 
reverse current reaches its maximum value. It is important to 
note that the turn-off process considered in this analysis 
does not imply ability of the thyristor to support forward 
voltage. 


Storage charge will depend upon forward current/ rate 
of change of current through zero, and junction temperature[3] , 
For a given, thyristor type, a range of values of storage 
charge will exist because of manufacturing tolerances. A 
cumulative charge distribution curve based on measurements 
on a sample quantity of a particular type of high power 
thyristor is shown in Fig. 3. 5 [3]'.' 
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The cumulative charge distribution curve can be approxi- 
mated by dividing the thyristors into a number of groups/ with 
all the thyristors in a group having equal storage charge. In 
Fig‘.'3.5/ group 1 consists of 10 thyristors (5 percent of the 
total) with = 300 micro C/^hile group 2 consists of 40 

thyristors (20 percent of the total) with Q =320 micro C/ 
etc". 

For convenience/ a second scale is introduced with 
origin as minimum value of 0j_^/ increment in 

v^ith respect to group 1 is designated as Q'/ the values 
corresponding to the groups of thyristors are = 0/ = 20 

micro C, etc". 

3. '2. 3 Combined Equivalent Circuits 

The combined equivalent circuit for the' supply system 
and thyristor valve is shown in Fig. 3. 6. The valve is consi- 
dered to be having N number of thyristors divided into n groups/ 
each thyristor having in parallel a grading network 
general/ the ith group of thyristors is represented by a 
switch Sj. / having in parallel a resistance series with 

a capacitance • 

3.2.4 sEimplification of the Equivalent Circuit; 

The equivalent circuit shown in Fig. 3. 6 is simplified 
for computof calculation and is shown in Fig. 3 .7, During the 
turn-off process the valve is represented by a time varying 
resistance (R) and capacitance (C). 
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^ Tu2r n*~0 f £ Ov03rvolt3Cfss 

Fig, 3 .4 is considered as the reverse current characteris- 
Lic of the thyristors of group 1. The storage charge of thyristgrs 
Or igroup 1 is neutralised at time t = T^/ and the thyristors of 
Liiis group turn-off which is simulated by opening the switdh 
Or Fig, 3. 6. At this instant the current flows through 
the voltage grading network of group 1, 


The initial conditions for turning-off of group 1 are 
calculated, using the rate of change of current of Eq,( 3 . 3 ), 
as follows. 


I 


o 



(3,4) 



I dt = 


V 

o 

2K 


From (3.4) and (3,5), 


(3.5) 


lo=/(wgr7L^ (3.6) 

and 

Tjl = /r2Q/(V^/L^) (3.7) 

After the switch is opened, the voltage across 
group 1 increases from zero as the current is diverted into its 
voltage grading network, while it continues to flow through the 
remaining switches ~ S^. After some time T^, the total 
current time integral will neutralise the storage charge of 
group 2, and S 2 is opened. Voltage will now start to build up 
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on the thyristors of group 2, as the current is diverted into 
its grading circuit. 

This process of opening of switches continues until 
all the thyristors are turned-off. As the group 1 turns-off 
first/ it experiences the highest voltage/ v^hile the group n 
experiences the lowest voltage. Starting with the turn-off 
of group 1/ a fast rising recovery voltage is goir.^rr.t -C , and 
is_ superimposed on the overvoltages caused due to non-simult- 
aneous turn-off of thyristors. 


3.3.1 System State Equations; 

The state equations are derived for Fig. 3. 7 using the 
algorithm given in Appendix A, These equations arewritten in 
the form 

px=Fx+Gu ( 3 . 8 ) 


where x is the state variable vector with 11 state variables 
given by 


[x] = [v,/ V / 


V. 


6 ' ' 7 ' '' 9 ' ^ 11 ^ 


^ 13 * ^ 15 ‘ 


18 ' 22 ' 3 


^ 19 ^ 

(3.9) 


F is a 11x11 matrix and G is a 11x2 matrix which are given 
in Appendix E.' u is a vector given by 

3.3.2 Computer Program ; 

The state equations are solved numerically using Runge- 
Kutta fourth order method. Trapezoidal rule ‘ 

... .. & 
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3.8 Flew chart or tlv.. coiaFCtor ;rocirrja 




integrate the current numerically for obtaining storage 
charge,' The flo'w chart of the computer program is shown in 
Fig. 3.8'.' 

An Exa mple 

3.4,1 Description of the Problems 

A valve vj'hich contains 200 thyristors and operates at 
a peak repetitive voltage of 200 k'V [ 3 ] is considered to 
show the effect of certain critical parameters on turn-off 
overvoltages with the following supply system and valve 
parameters'. 


Supply system: 


Valve: 



AQ 


= 0.05 to 0.3 H 
-2 nF 
= 2 nF 

= 1 m 

- 100 KV 

= 1 to 10 K-ohm 
= 10 to 500 nF 
= 5 to 30 ohms 
= 1 to 20 ' micro-F* 

=300 micro~C. 

= 100 micro-C. 


A valve consisting of 6 groups of thyristors is considered 
in this analysis. ■ Table 3.1 shows the number of thyristors 
in each group and their corresponding storage charges. 
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Table 3 .1 


Group No. of thyristors in Storage charge 
each group ^^rr niicro 

for each group 


1 

10 

190 

2 

30 

210 

3 

70 

230 

4 

60 

255 

5 

20 

275 

6 

10 

290 


The initial capacitor voltages are assumed to be zero. 

The initial conditions for the inductor currents are 

isCTi) = -20 Amps 
15 (Ti) = -20 Amps 
iig(Ti)= -20 Amps 

where T^ = 30 microseconds'.' 

The integration step size used in Runge-Kutta subroutine 
is O.'l microsecond. 

3.4'.:2 Results: 

Fig. 3"."9 shows the current and voltage waveforms during 
turn-off. Equations(3 . 6 ) and(3 .7) arc used to calculate the 
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current/ 1^/ and the time at which group 1 turns-off, At 

the time T^/ the impedance of the valve increases as the 
switch opens / and the current is diverted into the parallel 
grading circuit. The initial voltage at time across the 
grading circuit of group 1 thyristors/ is 

same as the voltage across the entire valve as the ocher groups 
are still conducting. The voltaffe across the first group of 
th/ristors increases after time T^ as the capacitors of the 
grading circuit become charged, although the magnitude of the 
valve current decreases. 

As the groups' are turned-off successively, the magnitude 
of the current decreases, but the voltage across each group of 
thyristors which has already turned off, increases'.' At time T^, 
the last group of thyristors will be turned off, and the 
recovery voltage builds up across the valve,' 

Effect of network parameters; 

The voltage developed across the first group of thyristors 
to turn-off, is a function of the current furnished by the- supply 
network. Investigation of the current distribution in a 
representative supply network during turn-off period leads to 
the following conclusions. 

The di/dt during turn-off is relatively low, and 
current flows through the parallel anode inductance As 

the current flowing through anode resistance is negligible. 
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R^can be neglected from the equivalent circuit [s] . Accordingly^ 
EqV(3.8) is modified and used in the analysis. 

Although the valve capacitance affects the turn-on 
overvoltages [l]'/ its effect is negligible here as the over- 
voltages during turn-off are controlled by the low impedance 
damping and ■ grading circuits. 

The effect of varying transformer induc-'.ance, on 

turn-off overvoltages is shown in Fig. 3. 10. For a particular 
set of the values of the parameters of the netv/ork and valve ^ 
increasing in increases the peak valve voltage, while 
reducing the value of dv/dt. 

Fig. 3. 11 shows the effect of damping circuit capaci- 
tance, on turn-off overvoltages. Both the valve voltage 

and dv/dt are reduced vjith the increase of damping circuit 
capacitance. But the very high values of are not 

preferable, as the losses in this circuit are proportional 
to the value of capacitance'^ 

The effect of damping circuit resistance, on 

turn-off overvoltages is shown in Fig. 3. 12. The value of 
dv/dt increases with the increase in R^- By observing the 
voltage waveforms for R^ — 1000 ohms, 5000 ohms and 10,000 
ohms, it can be seen that there exists an optimum 'for a 
particular value of at v;hich the valve voltage is minimum. 



65 


3.4.3 Discussion? 

The various results obtained are compared with that of 
ref, [ 3 ], and it is found that there is no significant difference 
whicn validates the program developed. In ref .[ 3 ] / anf IBM' 
supplied ECAP program is used for netv/ork analysis. Also the 
switching of individual thyristor group is represented by a 
step change in the impedance from a very small value to a 
large value. Hov/over, in this program^ i^'eal switch is used 
and the entire valve is represented by a variable impedance of 
R and C in series. This simplifies the equivalent circuit. 

The valve voltage calculation for a particular set of parameters 
is taken 19 seconds C.P.U. time on DEC-1090 computer, 

3 '. 5 Calculation of Recov ery Overvoltage Neglecting -Turn-o ff ' 
Process 

In the case of thyristor valves, because of the un- 
neutralized carriers (i.e. storage charge) present, a negative 
current flows through the valve during turn-off process, which 
delays the rise of reverse voltage across each valve. As it has 
been observed in the previous sections of this chapter, the 
recovery volt age across the valve is less dependent on 

storage charge distribution. Accordingly, the calculations 

\ 

in this section are carried out neglecting the turn-off 
process, i.e., all the thyristors of the valve are considered 
to be turned-off at the same Instant [ 4 ] , 
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Equation (3.8) is used for valve voltage calculations ^ 

The valve voltage variations v^ith time for different damping 

circuit resistances and capacitances are shown in Fig. S.lSo 

The results are similar to those obtained in the previous 

section considering the turn-off process. The frequency of 

oscillation and overshoot for particular values of end 

ate calculated which are useful for calculating the 

damping circuit losses considering an approximate second 

order system model used in the next chapter, fiith R^= 8K-ohm^ 
^d~ overshoot ( from Fig. 3.15) is approx. 25, 9 percents 

The expressions given by equations (3.14) and (3.15) [6j 

are used to calculate the damping factor of the oscillation and 

the frequency of oscillation respectively';' 

Mp = 100 exp(-Try )) (3.14) 

tp = 7r/[UJn ^Tl-Y^] (3.15) 

where Mp and tp are peak percentage overshoot and time 
to peak, respectively. 

For 25.9 percent overshoot, the calculated values 
damping factor and frejpiency of oscillation are 0';‘40 and 
4 '.5 KHz, respectively;' 

The results obtained in this section tally with the 
results given in ref, [4] without any significant difference. 



Time ( ji Sec) — — > 





voltage . ':■; 
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3.6 Conclus ions 

A .general computer program to calculate both recovery 
overvoltages and overvoltages due to non-simultaneous turn-off 
of thyristors in a.HVDC valve has been developed. The effect 
of varying the system and valve parameters on turn-off over- 
voltages has been demonstrated by considering a representative 
thyristor valve. 

/"although the effect of sv/itchyard capacitance, 
transformer capacitance, on turn-on overvoltages is significant, 
their effect on turn-off overvoltages is negligible. There is 
a significant increase in recovery overvoltage with the increase 
in transformer leakage inductance. 

Both the thyristor and valve voltages are limited by 
choosing proper damping and grading circuits. The overvoltages 
can be lilnited only by a grading circuit without any external 
damping circuit, but the values chosen for grading circuit 
resistance and capacitance are not feasible if we consider 
both turn-on overvoltages and losses into account. 

Consideration of turn-off proce5s in detail has no signi- 
ficant effect on recovery overvoltages across the valve, though 
its effect is considerable if we consider the voltage across each 
thyristor of a valve. A computer program has been developed 
for calculation of recovery oyrrvoltages neglecting the turn-off 
process'.' 
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CHAPTER 4 


ZJ- ® SENSIT I VITy CALCULATIONS USING 

ADJ OINT NEIVIORK APPROACH 


4.1 Introd uc tion 

In the design of any system, it is important to know 
the effect on the system performance due to variations of sorn^ 
system parameters. A method of calculating the sensitivities 
of the valve voltage with respect to damping circuit resis- . 
tance and canacitance in time domain using adjoint network 
approach is presented. The results for a representative 
valve are presented. This sensitivity information is useful 
in the design of valve damping circuit. 

4.2 Netv/ ork co nsid ered, a nd i ts Adj oint 

The equivalent circuit considered is the same as that 
presented in chapter 3 neglecting turn-off process. This 
network (N) is shown in Fig. 4,1 [4], The adjoint netvork N 
corresponding to the given network N is shown in Fig, 4 .2, 
which is constructed using. the algorithm given in ref.[5,ll]. 
The details of the algorithm are given in Appendix B, 

^ *3 S ensitivity Calculations 
4,3,1 System Equations: 

The state equations for the original network shown 
in Fig ,4,1 are given by Eq,(E.l) (Appendix E). The state 
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equations for the adjoint network are derived similarly and 
can be exoressed as 

p X = F X + G'u* ( 4 -,l) 

v/here x is the state vector given as 

X = Vg, V^3, V33, I3, ig, 

(4.2) 

F is a 11 X 11 matrix given by Eq.(E*3) of Appendix E. 

G* is a 11 X 1 matrix. The non-zero elements of G' are as 
follows S 

G*(l,l) = 1 “ [Ra/CR+R^)] 

G*(2,l) = ll - [Ra/(R+Ra)] (4.3) 

G‘(ll/1)= -RRg^/L^CR+Rg^) 

For studying the sensitivity of the valve voltage with respect 
to damping circuit parameters a current source (impulse 
function) is connected across R-C circuit of Fig*4»2V 

therefore u* is a vector c'iven by 

u‘ « [iQ(>^)] = (4.4) 

The initial state vectors for both original and adjoint 
network are zero .valued. 
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■Substituting the above expression for the u‘ in Eq.(4el)/ 
and integrating it between 0~ and 0*^ gives 

x(0’^) = -G» (4.5) 

Therefore^' solving the forced system Sq,’(4.1) is equivalent 
to solving an unforced system 

X = F X (4.6) 

with the initial state vector as -G*'.' 


4.3.2 Calculation of Sensitivities of the Valve voltage 
with respect to Damping Circuit Resistance and 
Capacitance s 


The perturbation observed in the valve voltage due to 
perturbations in the damping circuit parameters is given by 
the following equation (see Appendix B). 


AVo(tj) = s ^ [io (t) iR 

° ^ all O ^d ^d y'=t^-t 


s (v^ (r) v„ (t)dt] 


\T/ x'- / — - 

all C^o % ^d y=t^-t 


(4.7) 


The expressions for sensitivities are obtained from Eq.(4.7), 
wiiich are giv’^en by 




■■c>V^(t.)/aR, = s ^ [in ip 

° ^ all Rg o "d ’^d r =t£' 




(4.8) 


(4,9) 


all C^ o 


y =t£-t 
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Using Eqs,(4,8) and (4,9) ^ the expressions that are used for 
calcu3.ating the sensitivities in this analysis are developed 
ivhich are given by 






+ I iji^Cy') i^^(t) dt] 


(4.10) 


y = tj-t 


a vv 

ac^ 


t, 


■[ / V (t)dt + 2 v^rCV') V (t)dt 


13 'f ' "13 ■ “ ^ "15^' ' "15 

o 


Cm# j (>w 

V (t)dt] 

o 18 r = t,-t 


(4.11) 


Tho currents, voltages and the derivatives of the capacitor 
voltages, used in Eqs,(4.10) and (4,11) are obtained from the 
transient analyses of both original and adjoint netxvorks which 
are carried out numerically using Runge-Kutta fourth order 
method"," The integrals in Eqs.(4,10) and (4,11) are evaluated 
numerically using Trapezcidal rule [lO] , For example, the 
integral is approximated by the following siommation (with 
t^ = 5h) given by 
t 

I = [ i(y') i(t)]dt 

o ■ r = t£-t 

A ^ 

“lllT ~[i(5h)i(o)+2[i(4h^ i^h) + i(3h) i(2h) + i(2h) i^Eh) 

+ ^(h) i(4h)] + i(o) i(5h)] 


(4.12) 
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^ example 

4.4.1 Description of the Problems 

A valve which contains 200 thyristors and which 
operates at a peak repetitive voltage of 200 kv^ is considered 
to show the effect of damping circuit parameters cn sensitivi” 
ties with the folloidng supply system and valve parameters. 

Supply system; 

= 0.05 H 

C. = 2 nF 
t 

C = 2 nF 

y 

= 1 iWi 

d 

= 3 K-ohm 

OL 

V = 100 KV 
o 

Valve ; 

R , = 2 to 4 K-ohm 
d 

Cj = 50 to 100 nF 
^'d 

R =10 ohm 

.Kg 

C =1 microfarad 

g 

The step length used in Runge-Kutta subroutine is 0,5 micro- 
seconds. A C.P.U.- time of 42 seconds is taken on DSC-1090 
computer for the evaluation of the sensitivities, 

4 ■.'4.2 Results: 

The transient analysis of the original network (N) is 
carried out and the effect of damping circuit resistance (R^) 


\ 
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on valve voltage has been shown In Chapter 3 which is reproduced 

here in Pig',’ 4.3. The variation of the sensitivity of the 

valve voltage with respect to damping circuit resistance 

with time .for R , = 2 K~ohm and C = 50 nP is shown in 
■ ' o d 

Pig, 4, 4 7 It can be seen from Pig .4. 4/ that the sensitivity 
is negative from 0 to 160 micro secV, positive from 160 to 
440 micro sec,, and again negative from 440 to 600 micro sec. 

This variation of sensitivity can be validated from the valve 
voltage waveforms given in Fig ,4 .3 by comparison of the 
two curves corresponding to == 2K and R^ = 4K. 

The peak of the valve voltage for ^ K-ohm is 

occurring at 225 micro sec, (from Fig,4,3). The sensiti- 
vity corresponding to this time is positive (14,6 volts/ohm) , 

This indicates the increase in R^ will reduce the valve 
voltage at this time. The time to peak of the valve 
voltage is reduced with increase in r which confirms 
the results obtained in Chapter 3. 

The effect of damping circuit capacitance (C^) on 
valve voltage has been shown in Chapter 3 , which are 
reproduced here in Fig’*4', 5," The variation of the sensitivity 
of the valve voltage with respect to damping circuit 
capacitance (C^) with time for R^^ = 2 K-ohm and = 50 nP 
is shown in Fig.4-.'6t It can be seen from Fig ,4. 6, that the 
s:nsitivity is positive from o to 290 micro sec., and 
negative from 290 to S60O micro sec* 
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This variation of sensitivity can be validated from 

the valve voltage waveforms given in Fig. 4, 5, by comparision 

of the two curves corresponding to = 50 nF and = 100 nF'. 

The sensitivity corresponding to the peak valve voltage 

for = 50 nF±is positive (450 volts/nF) which indicates 

that the increase in reduces the valve voltage at this 

instant. Because of the negative sensitivity observed at 

the subsequent interval the time to peak of valve voltage 

is delayed with increase in C,. This is also validated from 

d 

the results obtained in Chapter 3'. 

Effect of and on Sensitivity: 

The variation of sensitivity of the valve voltage 

with respect to damping circuit resistance with time for 

different values of is shown in Fig ,4, 7. The variations 

in sensitivity of the valve voltage with respect to R^ with 

time for different values of are shown in Fig. 4. 8.' 

The sensitivity curves obtained for different values of R^ 

are similar in shape but the magnitudes are reduced as R^ 

is increasedi' It can be concluded that over a wide range 

of variation in R # the effect of increase in R, is to 

d ^ 

reduce the time to pesik of the overvoltage. The behaviour 
of 3 Vq/9r^ for different values of is different. The 
magnitude of the sensitivity is slightly increased in the 
first interval (when the sensitivity is negative) and 
reduced ih the siibsequent interval. The transition from 
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negative values to the posit! ve is also delayed and hence the 

peak overshoot is reduced with increase in from the initial 

value of 2 K-ohm, However# from the Fig*4«7# it can be observed 

that the peak overshoot is not necessarily reduced as 

Rj is increased •' For a given value of C,# there is an 

d 

optimum value of for which the peak overshoot is a minimum. 

The variation of sensitivity of the valve voltage v/ith 
respect to damping circuit capacitance for different c^ is 
shown in Fig.4,'9, The variations of sensitivity of the valve 
voltage with respect to with time for different Rd are 
shown in Fig74'.T0V The sensitivity curves obtained for 
different values of are similar in shape but the 
transition from positive to negative is delayed Vi^ith 
increase in C^. The magnitude of the sensitivity is also 
reduced with increase in It can be concluded as 

increases the valve voltage peak occurs at later instant 
and the effect on the peak voltage is diminished. The 
sensitivity curves for different R^ are similar in shape, 
but the increase in R^ reduces the magnitude of sensitivit;^, 

474 ',3 Calculation of Sensitivities of the Valve Voltage 
with respect to Damping Circuit Resistance and 
Capacitance considering Thrn-off Process : 

A detailed model which is presented in Chapter 3 
for considering the valve as a number of thyristor groups 
which turn-off sequential^ acxaording to their storage charge 
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distribution ^is incorporated in the sensitivity analysis.' 
Effectively the turn-off process is incorporated by 
increasing the values of grading circuit resistance (R) 
end capacitance (c) in steps^ till the last group of 
thyristors turns— off. Thus the system can be still 
represented by the equivalent circuit shown in Fig. 4.1 and 
the adjoint notvrork is similar to that shown in Fig.4"»2'» The 
only difference is the R and C are considered as variable 
elements depending on the turn-off process’;* 

The sensitivities of the valve voltage with 

respect to R, and C. for R . = 2 K-ohm and C, = 50 np 
d d d d 

are shown in Figs';'4.'ll and 4,12 respectively. The curves 
are qualitatively similar to those obtained in the previous 
section but the peak magnitudes are almost doubled. This 
indicates that the valve voltage variation with R^ or 
is more if the turn-off process is considered in detail. 

4,5 Conclusions 

A computer program has been developed to calculate 
the sensitivities of the valve voltage with respect to 
damping circuit resistance and capacitance using the 
adjoint network approach’.' The effect of damping circuit 
resistance and capacitance on sensitivities has been studied 
by taking a representative valve’;* 

From the sensitivity information it can be concluded 
that the effect of increase in advance the Instant 

of peak voltage whereas the effect of increase in is to 
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delay the valve voltage peakl The sensitivity is reduced 
at higher values of and C^. The sensitivities of the 
valve voltage V7ith respect to and have also been 
calculated considering the turn-off process in the valve. 
For a given value of R^ and it is observed that the 
S' nsitivity is increased compared to the case neglecting 


turn-off process. 
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CHAPTER 5 

CALCULATIO N OF D.AMPING CIRCUIT LOSSES 

5,1 In t roduc t 1 on 

The power loss Ip damping circuits is a significant 
component of the total losses in a convertor station. Also 
the amount of heat to be dissipated from the damper . raclshvr'^ 
affects their size, cooling requirement, and cost. A method 
of calculating the damping circuit losses has been developed 
representing the exact equivalent circuit of the convertor 
system during each interval of conduction over a crycle. 
Computer program is developed and the effect of bridge opera- 
tion (commutating and delay angles) and damping circuit 
parameters on losses is analysed using a numerical example. 

The results are comp'.j.'ed (1) with the results obtained using 
the approximate formulas given by Ainsworth [7], (2) and with 
the results obtained using an approximate second order model 
for simulating vailve voltage oscillations [4] , 

5'.2 Ideal Valve Volt age , Wave form . 

The three phase bridge circuit (shown in Fig. 3.1) is 
considered. In the normal operation of the bridge circuit, 
two or three valves conduct alternately . Considering the 
firing of valve 5 at the instant a, neglecting transients. 
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the voltage appearing across the valve 3 over a cycle is 

given in Table 5.1^ with, the supply voltages given by Eq.{5ol)[4] . 

v^ = 1^ sin(o^t + 90°) 

v^ = El sin(Cht - 30°) (5.1) 

V = /“ E_ sin(CX?t ■“ 150°), 

C 1 3 L 

where E_ is the line-to-line r..m.s. voltage, 

Ij 

The valve voltage x\?avefonn for a -^= 15° is shown in 
Eig *5.1. 

5V3 Ainsworth* ^_Formulas for Damping Circuit Loss Calculati^n [ 7] 

It is assumed that the current in the damper resistor is 
that obtained by applying the ideal valve voltage shown in 
Fig. 5.1 directly to the R-C circuit, neglecting the transient 
voltage drop in the commutating reactance. One cycle of the 
waveform of the ideal valve voltage consists of eight arcs 
of fundamental frequency sine waves, . :separated by vertical 
discontinuities or jirnips — one such arc (for the period of 
conduction of the valve) being of zero amplitude. The curient 
resulting from this voltage wave in each of the eight intervals 
corresponding to the separate sinusoidal arcs has two components: 
(a) a transient component i^ resulting from the preceding 
voltage jump and (b } a sinusoidal component i^ resulting from 
the sinusoidal arc of voltage* 
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Table 5.1; Definition 

of the 

Valve Voltage Trjaveform 

Sec- 

tion 

Time Interval 
(in radians) 

Valves 

conduc 

ting 

Theoretical vo.ltage across 
valve 3 , v (60 t) 

1 

a < < a + /-c 

3,4,5 

Vj^ = 0 

2 


4,5 

^2 ~ 

3 

j -KC <CA?t < j 

4,5,,6 

V_ E_sin(C*2t4 ? ) 

3 \ 2. Li 6 

4 

j + < ^ + a 

5,6 

V 4 = V 2 sin<;t>t 

5 

~ + a <6^t < 4a+/^ 

5,6,1 


6 

™+a+^<cJt < 71 4a 

6,1 

Vg = Y 2 E^sin( 6 u)t- |) 

7 

7r4a < 0 >t < 

6 , 1,2 

Vrj ~ E^sinCi^t- ~) 

8 

1 47T 

7 T 4 a 4 ^<^t < “ + a 

1,2 

Vg = E^sin(a;t- |) 

9 

r» 14. ^ 471,^ , . . 

^4- CC <(^t < “'-4a4yMv 

1,2,3 

Vg = 0 

10 

|2 -Hi+^<CJt < ^ + a 

2,3 

^10 = ° 

11 

^*—4- a ■t^^t 4 4a4'yf4.. 

2,3,4 

Vii = 0 

12 

™ 4- a4«l.<4Jt < 27T4a 

3 ^ -4 

3,4 

^12 




Fig. S.lTheoretical voltage acres? rectifier valve, showing voltage jumps. Drawn for a 
««« 15": Vji a is ignition voltage jump; Fj,=-\Si%jin S is extinction 

voltage jump. The smaller jumps (F^i/2 and occur at ignition and extinction of 

other valves. 
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The sinusoidal component of current is an arc of 
fundamental frequency leading the corresponding voltage ■ 
wave by 90 ^ since at this frequency the reactance of the 
capacitor predominates. The instantaneous power consum.ed by 
the resistor is given by 

p = R i^ = R(if + 2i i + i^) Watts (5‘.2) 

t t s s 

In the above equation/ the first term is the high 
freque'ncy loss/ and the sum of the second and third terms is 
arbitrarily called the low frequency loss/ although the second 
term depends on both components of current. It is obvioaas 
that the first and third terms are always positive but that 
the second term may have either sign. The expression must be 
integrated with respect to time over the eight parts of one 
cycle; and the result must then be multiplied by fundamental 
frequency in order to find the energy dissipated per second, 
which is the pov/er loss. 


Ainsworth has derived the loss formulae which are 
as follows: 


where 




p. , = 1.75 f.vf (C4C )(sln a + sin 6) 
hf 1 LL s 


(5.3) 


(5.4) 


2 2 2 

p = 2Trf., v__ C R[2V46+0.875(sin2a + 3 sin26-2u) 

Xjjl J- LiLj 

+ O'.433(cos2a + cos25)] 


(5.5) 
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= fundamental frequency ^ Hz 
Vtt “ r.rn.s. llne-to-llne secondary voltage> volts 
C = damper capacitance, farads 
Cg = stray capacitance across the valve, farads 
R = damper resistance, olirne 

a = ignition delay angle ; 6 = extinction delay angle, 
u = overlap angle, radians 

= power loss in damper for one valve, watts 
= high frequency component of P^ 

Pj^^ = low frequency component of P^. 

High Frequency Loss; 

The voltage jump V . applied across an R-C circuit 

2 

causes a loss of energy of CVj/2 joules in the resistor. The 
voltage across a valve has eight jumps per cycle, the largest 
of I'vrhich are the iqnj.tion jump V . . = /2V^ sina and the 

J 1 j-jIj 

extinction jump V._ = V sin6. In addition, there sus 

J 0 Jui-i 

are six smaller jumps due to comrnutations in other valves; 
three of these are equal to dnd three are equal to 

V . /2V The energy loss per cycle i^, therefore, 

J® 

/ 

W, - = 0.5C S = 0.5C(\r2 V__)^ (sin^a + sin^6)[l+3 (0V5)^] joules 

inx; j Lij-i 

(5.6) 

which multiplied by f^/ gives the high frequency power loss 
as given in Eq*(5','4). 
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In the total power loss, high frequency lo^ constitutes 
a major part and as an approximation the lov/ frequency loss may 
be ignored! The high frequency loss is independent of the 
resistance and is directly proportional to the fundamental 
frequency, to the capacitance, and to the sura of the squares 
of the voltage jumps. 

^ '• Calculation of Da mping Ci r cuit L osses usi ng a seco nd 
Orde r M o del f or Vol tage O s cillations [4] . 

After each commutation, the exact valve voltage consists 
of a damped oeillation vjhich is superimposed on the sinewave. 

The frequency and amplitude of the superimposed damped 
oscillation are calculated as explained in 3-. 8 of recovery- 
voltage calculation in Chapter 3. In this analysis, it is 
assumed that the same type of damped oscillation is superimposed 
on the ideal valve voltage at each commutation. 

The equivalent circuit used for power loss calculation 
is shown in Fig. 5. 2 uhere the valve voltage is generated by 
svjitching one of the generators (v^ - of Tab'le 5.1) to a 
linear seconrl order system which superimposes a damped 
oscillation on the sinev/ave. The output voltag-e of the 
second order system is applied to the damping circuit. 

The loss in the damping resistance, R^, is given by 

T 

p = (rVT) / i^ (t)dt (5.7) 

o '^d 
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wh'^re • instantaneous current through 

R^ : damping resistance 

T ; power frequency period. 


The transfer function of the second order system, is of the 
form 

Y(s)/u(s) s +OJ^) (5.8) 

where 

: damping factor 

Cl).^ i undamped natural frequency. 

Choosing the state variables as x^/ (x^ = Y, X 2 = dY/dtj,)^ 

where y is the output of the second order system), and the 

capacitance voltage v^ , the system eauations are developed [&] 

d 

and are given by 


p^l ’ ■ 


1 PX2 


PV^ 


L_ 1^-' 



1 0 1 


r — 

^1 


_ 

0 

2^0^ 0 


^2 



0 --1/C,R, 

d d_ 




0 


^Rd “ 


[U] 


(5.9) 


These e^^uations are solved numericali^using Runge-Kutta 

i '' 2 . 


fourth order method'.' After Icnowing the state variables ,, x 


and V at any time instant t, the current through the 


cl 


damping resistance (R^) can be evaluated'. Trapezoidal rule 
is used for performing ’the numerical integration required 
in Eq.'(5.7) of power loss calculation. 



97 


^ A. Anal ysis for D ampin g Circuit Loss Cal culation 
5.5.1 Detailed Equivalent Circuit Representations 

In the previous section, the ideal valve voltage is 
superimposed with a damped oscillation at each commutation. 

This is generated using a fixed second order model which is 
an approximation. Using the exact waveform of valve voltage 
is necessary to calculate the damping circuit losses more 
accurately.' This can be achieved by representing an exact 
equivalent circuit during each interval of conduction over a 
cycle and carrying out the transient analysis',' 

The equivalent circuits of the thyristor valve converter 
system, during each interval corresponding to sections 2 to 9 
of Table 5.1 are developed'.' The ec[uivalent circuit during 
the interval < (7T/3)+a (valves 4 and 5 conducting) 

is given in Fig. 5. 3.' In deriving this equivalent circuit, 
the conducting valves are represented by short circuits, and 
the non-conducting valves are repres'"nted by their corresponding 
gra.'-Ung circuits in parallel with their danping circuits. The 
anode damper R -L is neglected for all the valves except 
for valve 3 . 

The state variable equations are written (using the 
algorithm in Appendix A) to the equivalent circuits developed 
for sections 2 to 9 of Table 5.1 and are solved numerically 
using Runge-Kutta fourth order method. The final state vector 







99 


3t ths snd of S3ch intsjrvsl is tsJcon 3s tliG initisl st 3 t© 
vector to the subsequent interval. The current flowing 
through of valve_3 during each interval is evaluated. 
Trapezoidal rule [lo] is used for performing the numerical 
integration required in Sq,(5.7)"' 

5.5.2 Simplified Equivalent Circuit Representations 

The equivalent circuit of Fig. 5. 3 can be simplified 
by neglecting capacitances C and c and the anode danper 
(R^ - L^) of valve 3. The state equations for the simplified 
equivalent circuits are developed and solved numerically; 
and the power loss in the damping circuit of a valve is 
evaluated as explained earlier. The equivalent circuits 
for sections 2 to 9 (see Table 5.1) and the state equations 
for the first two intervals (sections 2 and 3 of Table 5.1) 
are given in Appendix C. 

5V5.3 Computer Program; 

The flow chart of the computer program is shown in 
Fig.5.4".‘ 

An Example 

5".'6'."1 Gy«tem Data: 

Supply System; = 0,05 H 


c^ = 2 nF 
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Fig'.'B z Flow chart of the computer program. 






10,1 


Valve : 



2 np 
1 niH 


= 3 K-ohm 

cl 

Et = 200/V2 Kv 


f = 60 HZ 


R, =2 K-ohm 
ci 

= 100 np 

Rg = 10 Ohm 

C =1 microfarad 
<3 

N = 200 


5 . 6V2 Results : 

Computed valve voltage, current in the danping circuit, 
and the damping circuit loss for a = ^0°, jU. = 15*^ are shown in 
Pigs. 5.5, 5,6 and 5.7 respectively, it can be seen from the 
curves that at each commutation, the valve voltage jump with 
a subsequent oscillation v/hich in turn produces a current 
oscillation, the magnitude of which is- relatively large 
compared to power frequency current. It carj. be seen from 
the power loss curve, loss increases rapidly at each commu- 
tation while varying slightly between commutacions , This 
confirms that the lo.sses are produced mainly by charging and 
discharging of when voltage jumps appear across it. 

Effect of Commutation Angle (f<.) and Delay Angle (a) on 
Damping Cii'cuit Loss: 

The effects ,of commutation and delay angles are 
shown in Pigs." 5.8 and 5. 9, In Fig. 5, 8 the solid lines 








Damping circuit loss (K.W,) 



cot (degrees) 

Fig. 5.7 Damping loss vs time 




Using the formulas given by ~Ainsworth 

Considering a second order system for oscillafior 
(50% overshoofyfrequcncy of oscH{ation=5KHz) 


= 75 


.o( = 90' 
oC = 60 



Commutation angle (Jli) 

Fig. 5.8 The effect of the commutation angle (| j) and 
delay angle on the power loss • 
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represent the calculations by Ainsworth's formulas/ whereas 
the dashed lines represent the calculations using a second 
order system for oscillations. It can be seen from Fig 5,8 that 
the increase in a increases losses/ because for high a the 
heights of the extinction and ignition voltage jumps which 
appear across the valve are more/ v/hich in turn increases the 
peaks of the current spikes. The increase of losses with^ 
is less significant compared to that of a. 

The effect of ^ on losses for two different values of 
(25° and 90°) for the three methods described is shown in 
Fig';'5'.'9. 

Effect of Damping Circuit Parameters on Damping Circuit Loss: 

The effects of and C^ on damping circuit losses are 
shown in Figs, 5,10 and. 5,11 respectively. It is observed 
that losses increases drastically with increasing C^/ whereas 
does not significantly affect the circuit losses. The 
capacitance C^ must therefore be kept to a minimum to reduce 
losses/ whereas the selection of the value of R^ depends mainly 
on the recovery voltage transient, 

5.6I3 Discussion: ' 

Three methods have been studied for calculating the 
damping circuit losses. In the case of Ainsworth's formulae 
an ideal valve voltage waveform is considered for loss 
calculation neglecting the valve voltage oscillations during 







0. 

50 

TOO 

Cd (nF) 
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200 

Fig . 5 .ir 

Effect 

of Cd on 

damping 

circuit loss* 
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commut ationV In the second method, an improvement is made by 
assuming a fixed second order model for simulating high frequency 
valve voltage oscillations, that are imposed on the ideal valve 
voltage. A detailed analysis is carried out in the third method, 
in which the exact equivalent circuits are represented during 
different intervals to evaluate the exact valve voltage waveform 
Xf-rhich can be used for calculating the losses more accurately. 

It is observed that Ainsv/orth's formulae give inaccurate 

results at lower values of a both qualitatively and quantitac- 

ivelyV From Fig. 5, 10, it can be s^'en that increase in gives 

in 

rise to rapid increase^losses at lower values of which is not 
true? and also the losses arc almost constant with increase 
in at higher va3ues of a, although the exact losses are 
decreasing. At higher values of a the formulae give reasonably - 
accurate results. 

In the method considering a second order model, 
reasonably accurate results can be obtained (though the 
deviation is considerable at higher values of Ct)^provided the 
overshoot is chosen correctly. In Fig. 5, 9, the 'curves with 
two different overshoots, 50 percent (v/hich is assumed), 

25.9 percent (which is obtained from the transient analysis) 

! 

show that the results can be improved by considering the over~ [ 

shoot obtained from transient analysis. It is to be noted that [ 

i; 

the losses increase significantly with increase in the over- | 

r 

!; 

shoot, whereas the effect of frequency of oscillation on losses j 
is not significant [4]‘, ; 
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A Gomparision In terms of accuracy and computational 
time for different methods is shown in Table 5'.2V 

It can be observed from Table 5 , 2 , the step length 
required in the detailed analysis is half the step length 
rer-uired in the second order system model. The detailed 
analysis, of course, requires much more computational effort 
than the first tvro methods and is used here mainly to investige.te 
the accuracy of the earlier methods. 

^ Cpnc fusions 

Three methods of calculating the damping circuit losses 
are studied.' An exact calculation of the circuit loss is 
carried out by using detailed e<fuivalent circuits during each 
section of the valve voltage waveform. 

The effects of varying the bridge operation parameters 
(ct andyW. ) and the damping circuit parameters ^nd C^) 

analysed for a representative thyristor valve. The damping 
circuit losses are increased by but less affected by 
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Table 5V2: Comparison of various methods used for 
damping circuit loss calculation. 


Method 


Damping 
circuit 
loss for 


Step length 
used 

( seconds) 


C«P.U.time 
on DEC- 10 90 
computer 
( seconds ) 


15 


Using Ainsworth’s 404.659 ; i- 0.09 

formulas 


Using second order 
model for oscillations 
(overshoot 25,9 per- 
cent, frequency of 
oscillation 4.5 KHz) 


423.6 0.15625x10“^ 


5.07 


Using second order 
model for oscillations , 

(overshoot 50 percent , 507.11 0.15625x10“^ 5.07 

frequency of 

oscillation 5 KKz) 


Using detailed analysis- 

simplified ecfuivalent 416,03 9 0',78125xl0 3 9,04 

circuit representation 


Using detailed analysis-^ 

detailed equivalent 391,442 0,78125x10 

circuit representation 


55,00 



115 


CHAPTER 6 

VALVE DAMPING CIR CUIT DESIGN FOR HVDC SYSTEM S 

6 VI I ntroduction 

The electrical design of HVDC valve depends upon the 
overvoltages generated in the system and those generated 
xvithin the valve. These overvoltages will depend upon the 
valve pararneters. In the case of thyristor valves^ the voltag 
grading circuity the firing method/ and the mechanical 
design are the most important factors. The overvoltages 
developed during the turn— on and turn-off periods are .among 
the most severe internal overvoltages. The valve damping 
circuit which consists of a resistance and a capacitance in 
series is used to limit the recovery voltage during turn-off. 

One of the important factors to be considered in the 

design of valve damping circuit is the damping circuit loss 

that takes pi ace because of the voltage applied across the 

the 

damping circuit', when/valve is not conducting. Proper 
selection of damping resistance and capacitance limit feHe 
recovery voltage to an economic value without undue 
increase in the losses of the damping circuit. Calculation 
of turn-on and turn-off overvoltages has been presented in 
Chapters 2 and 3 respectively. Calculation of recovery 
voltage and damping circuit losses has been described in the 
Chapters 3 and 5 respectively. In this chapter,^ a method has 
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been proposed to optimize the damping circuit with the objective 
of limiting both recovery voltage transient and damping circuit 
loss . 

6.2 Erru i vale nt Cir cu i^ 

The equivalent circuit used in this analysis- is the 
same as the equivalent circuit used in the case of turn-off 

i 

overvoltage calcula -tions (which is given in Fig. E.l of 
Appendix E) except' the following changes; 

(1) The source voltages are taken to be sinusoidal 
voltages . 

(2) The turn-off process -is neglected. That is, the 

values of R- and C in are set equal to R_ 

f' (jr 

and Cq respectively. 

^*3 A Method for Parameter Optimization 
Given a system described by 


i = A(p) X 

where p = (p^^, p^) 

p^,i=l,2y . . . /n are the system parameters',’ 

the problem of optimization of the parameters p, for the 

above system is defined by the objective of minimizing an 

index of performance^, given by 

0^= x"^ Q X dt (6.2) 

o 

where Q is a positive semidefinite matrix. 
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The criterion for selecting a performance index (PI) 
is dependent on the desired objective to be achieved and is 
accordingly expressed as a function of the design verifies « 

It can be shown that for any asymptotically stable system^^ 
the PI can be expressed as given 

^ = ^o P "^o 

where p, a symmetric positive semide finite matrix for any 
stable system# ds given by the matrix Liapunov equation 

a'^p + pa = .-tQ (^?4) 

x^ is the initial state vector. For a given set of initial 
values of parameters pw...#p # the matrix A is determined 
and the matrix Liapunov equation is solved. The matrix P 
is used to evaluate the PI. By searching through the parameter 
space, it is possible to determine the optimal values of 
pararaeters that minimize the PI. 

A solution procedure of the matrix Liapunov equation 
(Eq. 6.4) is given in Appendix D. 

6.4 Stcat e e quations 

The state equations for the equivalent circuit 
considered, is given in Appendix E. These are of the form 
px = Fx + Gu (^^5) 


where u is a vector given by 
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where 


u = [v, v_] 


V^ = V, ~ V 

1 b a 


-2 "b ~ 


( 6 . 6 ) 


(6«7) 


In order to use the method described in the previous 
section^ the Eq,{6„5) should be transformed to the form of 
Eq*(6.1)';' Since the vector u consists of sinusoids] voltage’ 
direct transformation will make the system critically stable 
and the technique cannot be used.. To avoid this difficulty- 
damped sinusoidal voltarfes are assumed. The damping term is 
chosen arbitrarily. A secozid order system of the following 
form 



7 n dt n 


0 


gives the solutio'n 

y(t) = emp(‘^C*^t) sinCdt 


with the following initial conditions 


y(o) = 0 
y (o) ~ (X> 


Taking the state variables 
respectively/ the expressions for v^^ 
and x ^2 given by 


13 

and 


'.■<,) ar5 -?(t) 
in terms of x. 
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V, 




2 V3 ' CD' ^12 " ^13^ 


''2 =\/f 


OJ 12 ^^13 


Eq7 (6.5) is transformed to an unforced system of the form 

p Xt^ = 

where the state vector x^ is given as 

“ [^6' ''7' ''9' ''ll' ''13' "15' "18' "22' ^^3 ' S' ^9’ 

lT 

^ 12 ' '"l 3 ^ 


Fj^ is a 13 X 13 matrix with the following non-zero elements 


Fj^(i,j) = F(i,j) for i^ j- = 1,2,,.. ^11 
Fi(9'12) = [/I 1 /3)]/3L^ 

Fj_( 9,13) = [^ E^(l/Cj)]/3L^ 

F3_(10,12) =[^ -V3)]/3L^ 

F^(10,13) = -^1 E^(l/Cj)/3L^ 

F 2 ^( 12 , 13 ) = 1 
F 3 _( 13 , 12 ) = - 0 )^ 

F^(13,13) = 
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6’,'5 Selection of the Performance index and Initial Vector 


The performance index chosen for limiting both recovery- 
voltage transient and damping circuit loss is given by 


or 


o 

OL 


2 


^ T 

=i /• X Q X dt 

(6.13) 

o 


T 

dt = / x Q, X dt 

o 

(6.14) 

= / X X dt 

(6.15) 


o a o 

e is the transient voltage (appears across the valve )which 
is given by 

e = X (6,16) 

■whore 


C = f — — 2— 

^1 ^R^+R R^+R. 


R 


0 0 0 0 0 




^ ° Al- ° °] 


'•G a "'G a 

(6.17) 

e^ is the steady star.e voltage across the valve 'which 
by 

(6.18) 


®a = =2^ 


v;here 


c, = [00000000000 ^2E O] 


(6.19) 


i„ is the current through the damping circuit which is given 

R -j 

d 

by 


^R 


= D X 


( 6 . 20 ) 


d 


v/here 

D = [l/R(a 1/R^ 0 0 0 0 ~1/R^ 0 0 0 0 0 O] (6.21) 

K is weighting factor given to the po'werloss component. 
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The above performance index consists of integrals 
that are to be evaluated from zero to infinity. Although, 
the damping circuit loss is be considered only during a 
finite period, it is assumed here that minimizing the Pli ^2 
is equivalent to the minimization of the actual povmr losses'. 

The weighting matrix C is given by 


0 = 

% 

K 

(6.22) 



t" 


Qi = 

(Cl - 

■ ^2^ ^^1 " ^ 2 ^ 

(6.23) 

ii 

a 


D R^) 

(6.24) 


The initial state vector considered, corresponds to the instant 
a = 90°, and is given by 

Xq= [0 00000000000 ^ 6 , 25 ) 

An Exa mple 

6.6'.l Description of the Pro'JDlems 

A valve which contains 200 thyristors and operates at 
a peak repetitive voltage of 200 KV, is considered v/ith the 
following supply system and valve parameters,'. 

= 0.05 H 
= 2 nF 
C =2 nF 

y 

R =3 K-ohm 
a 


Supply system: 
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L = 1 mH 

3 . 

= 200/^2 KV 

ij 

f = 0.05 
f = 60 HZ 

Valve: = 2 K*-ohm to 8 K-ohm ? C, = 50 nF - 200 nF 

a d 

R =10 ohms 

g 

C =1 microfarad^ 

6 '. 6 . 2 Computed Results and Discussion; 

Fig. 6.1 shows the variation of Pi with for 
== 100 nFV lAihen K = 0, the optimum R^ (at which PI is 
minimum) is obtained as 6.0 K-ohm, 

As verification to the optimum R^ obtained for 
= 100 nF/ the recovery voltage transient for different R^ 
is computed from the transient anajysis of the e-uivalent 
circuit (v/hich is explained in 3,5), The effect of R^ on the 
recovery voltage transient is shown in Fig. 6 , 2, It can be- 
seen from Rig. 6 . 2 that the recovery voltage transient is 
mnlmura for = 6 K-ohm comparod to that obtained with other 
values of R^V Though the peak overvoltage is less for R^=4K-ohm/ 
the deviation of the transient voltage from the steady state 
voltage is more and persists longer (whereas for R^ = 8 K-ohm/ 
the deviation of the transient voltage from the steady state 
voltage is less but the overvoltage peak is more compared to 




Fig. 6.1 Variation of PI with for C<i ^lOOnF * 
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= 6 K-ohm case),' Thus It is observed for = 100 nP, 

= 6 K-ohm gives the minimum recovery voltage transient'. 

For different values of the optimum is evaluated 

and the Pi values corresponding to optimum R^ are given in 

Table 6,1. it can be seen from the table (with K = 0) that the 

optimum value of is around 100 nF. Therefore ^ for limiting 

the recovery voltage transient the optimum damping circuit 

parameters are R. = 6 K-ohm, C, = 100 nF. 

d d 


Table O','! PI Values corresponding to optimum R^ 
for different C^. 


^d 

Optimum R^ 

PI 

50 nF 

6,5 K-ohm 

O'. 574 543x10^ 

100 nF 

6,0 K-ohm 

0.533494x10^ 

150 nF 

6 .0 K-ohm 

0.555074x10^ 

200 nF 

6,0 K-ohm 

0.594936x10^ 



3 


dptimj._ dampina cir cuit par a meters with inclusion 

The effect of the inclusion of the power loss in the 
PI is shown in Table 6.2 which gives the optimum values of 
and for different weighting factors K. 

Table 6,2; Optimum damping circuit parameters 
for different weighting factors(K). 


Weight age given Optimum Optimum 

to power loss 


1 

5*5 

k-ohm 

100 

nP 

10 

6,0 

K-ohm 

50 

nP 

50 

4.5 

K-ohm 

50 

nP 


^ ^ Con clusi ons 

The design of the damping circuit requires the deter- 
mination of the capacitance and resistance v;hich reduce the ' ! 
recovery voltage transient without undue increase in losses 
in the damping circuit. An objective function is derined which 
• takes into account both the recovery voltage transient and 

the danping circuit loss* The optimum damping circuit parameter; 
are evaluated for a representative valve using the approach of ,, 
minimization of the objective function. 
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CH/lPTER 7 
CONCL U SIONS 


7Vl Summary 

The electrical design of HVDC valves depends upon the 
overvoltages generated in the system and those generated i;ithin 
the valve. The overvolt ages / developed during the tum-on and 
turn-off periods/ because of the non- simultaneous operation of 
the thyristors in a valve/ are among the most severe ’ over- 
voltages'! The grading and damping circuits are used to limit 
these overvoltages instead of derating the valve! The design 
of damping circuit also depends upon the losses that take 
place in it! Proper selection of damping resistance and 
capacitance limit the overvoltages to an economic val<^ without 
undue increase in the losses of the damping circuit. 

In Chapters 2 and 3 of this thesis / an attempt has been 
made to calculate the turn-on and turn-off ;overvoltages in a 
HVDC valve! Computer programs are developed and the effect of 
various parameters on these overvoltages has been studied. The 
approach adopted for simulation of a valve during turn-on and 
turn-off periods is along the lines given in refs, [l] and [s], 
but the equivalent circuits are simulated by developing the state 
eauations on the basis of topological considerations/ instead 
of using the ECAP program for simulation as in refs, [l] and 


[ 3 ] 
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In connection with damping circuit design^ a computer 
program is developed to calculate the sensitivities of the valve 
voltage with respect to the damping circuit resistance and 
capacitance in time domain using adjoint network approach using 
the algorithm given in refs, [5, ll] / in Chapter 4. 

In Chapter 5, a detailed method has been developed to 
calculate the damping circuit losses by; .representing exact 
equivalent circuit of the converter system during each interval 
of conduction over a cycle. The results obtained are comparable 
(1) with the r^'- suits obtained using a second order system for 
valve voltage oscillations [4]/ and (2) with the results 
obtained using the approximate formulas given by Ainsworth', 

In Chapter 6^ an objec.tive function, which considers 
both recovery voltage transient and damping circuit loss, is 
defined and optimal damping c.ircuit parameters are obtained by 
minimizing the objective function. 

S cope for Fu r ther Work 

The calculations of turn-on and turn-off overvoltages 
can be greatly improved by considering the number of thyristors 
which turn-on and turn-off non- simultaneously as a random 
variable and then evaluating the expected value and variance 
of the overvoltages, instead of considering the groups of 
thyristors which turn-on and turn-off at a specific time. 



127 


The sensitivities of the valve voltage with respect to 
damping circuit resistance and capacitance can be used for 
optimizing damping circuit parameters by defining an error 
function . The gradient of the error function can be ascertained 
and the parameter values can be altered according to some 
gradient technique to efficiently reduce the error| 

The technique presented in Chapter 6, for optimization 
of damping circuit parameters is a first step in the application 
of mathematical optimization methods to damping circuit design, 

( The assumptions used in applying this technique have to be 
examined critically^ particularly regarding the choice of a 
Xz suitable performance index. This can lead to improved method 
selection of damping circuit parameters. 
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APP-RMDIX A 


FORMULATION OP ST'\TE EQUATIONS FOR LINEAR 


networks [ 12 ] 


Consider a network containing only two terminal 
components eacli of tvhich is either a voltage (or current) 
source or is a passive element (resistor, or inductor, or 


capacitor) . The voltage current 

relationships of 

the -passive 

elements are 



v(t) = R i(t) 

for resistors 

(A.l) 

i(t) = C(dv(t) /dt) 

for capacitor5 

(A. 2) 

v(t) = = L(di(t)/dt) 

for inductor5 

(A.3) 


where v(t)_ and i(t) represent the voltage and current variables 
of the elements'." 

If the edges in the netvrork graph corresponding to 
the voltage souarces do not form circuits and the edges corres- 
ponding to current sources do not form cutsets, then the 
network can be characterized by a set of linear differential 
eauations vbich can be derived as follows. Let the voltage 
and current variables associated with the edges of an 
arbitrary tree T be designated by the vectors v^ and i^. If 
the vectors Xj^ = ±^) and represent. 
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respectively y, the primary and secondary variables corresponding 
to the tree T and the cotree T*/ it is possible, by proper 
partitioning of these vectors, to write the terrdnal equations 
for the components in the form 

= specified functions of time (voltage and current 


sources') 


d 


■dt%i = 

and 


(A.4) 




^3 


(d/dt) 


1 


0 


D_ 


L 




s2 




S3 


(A. 5) 


Although it is not evident from the development so 
far, it can be seen later that the vector ^ep^^ssents the 
state vector of the system. The matrices Dj^(i = 1,2,3) 
in Eqs.(A,4) and (A.5) are all diagonal with positive diagonal 
entries, and corresponds to the resistance eleraents',' 

The vectors of primary variables i = 0,1, 2, 3 
contain exactly one variable for each edge of the graph and 
are identified with the voltage and current variables of the 
various elements’.The vectors of secondary'^'^^i^^l®® ^ ~ 

2,3, likevjise contain exactly one variable for each edge of the 
graph and represent the complimentary variables of the respective 
primary vectors’.' The circuit end cutset equations of the system , 
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graph for any tree T can be written in the form 








= 0 (a. 6) 

where each one of the submatrices S6 . . has the form 

J 



where A . . is the fundamental cutset matrix and B. . is the 
fundamental circuit matrix, 

A maximally selected tree and cotree by definition 
include a maximum number of edges correspcnding to capacitors 
in T and a maximum nuimber of edges corresponding to inductars 
in T*'. It can be shovrni/ for such a tree, certain of the sub~t 
matrices in Eq.(A,6) are identically zero and that the vector 
of primary variables"^^ identified with the tree and cotree 
represents the state vector of the network’.' 

A systematic procedure for identifying a maximally, 
selected tree in the network graph G is based on a sequence 
of subgraphs of the network graph selected as follows* 
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Consi''er a subgraph of G consisting of all the edges 
corresponding to voltage sources’.' Select a tree in G^‘;‘ 
By a hypothesis o f the theorem^ the voltage sources form no 
circuits. Consequently^ in G^ the coi'npllTaent of 
contains no edges’. 

Consider a subgraph G2 of G consisting of all the edges 
in Gj and the edges corresponding to capacitors. Select 
a tree in G2 containing in G2, the cov»5?lixT:ent of 

T2 contains at most edges corresponding to capacitors, 
i.e. the edges associated with the voltage variables 

Consider a subgraph G^ consisting of all the edges in 
G2 and the edges corresponding to the resistors, i.e, 
the edges corresponding to the variables in , 

Select a tree T-, in G_ containing T . in G_, the 
complement of contains at most edges corresponding 
to capacitors and resistors, i.e., the edges corres- 
ponding to the voltage variables in 
respectively',’ 

Consider a subgraph G^ of G consisting of all the edges 
in G^ and the edges corresponding to inductors,' Select 
a tree in G_^ v/hich includes . In G^, the compliment 
includes, at most edges corresponding to capacitors, 
resistors and inductors, i.e., edges corresponding to 
the voltage variables in’Xpg2''%3' ^'^'^'^1 


in T}^2 current variables in ZjJ 


p 2 
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5.' Consider the entire graph G. Select a tree = T which 

includes T^. Since, by hypothesis, the graph contains no 
cutsets of edges corresponding only to current sources, 

Tg - = T* and the maiciinally selected tree and cotree 

have been identified. 


From the second step in selecting the tree, it is 
evident that the fundamental circuits defined by span only 
those edges of the graph that are contained in T 2 » thus 
follows that the fundamental circuit equations introduce only- 

zero entries in ^22 ^73' ®22 ~ ^23 ~ Like- 

wise the circuits defined by the edges of span only those 
edqes of the graph contained in T^/ and the fundaiTiontal circuit 
equations can introduce only .zero entries in ^^ 2 * 

®32 “ 


From Step 4, it follows that the cutsets defined by the 
edges of corresponding to inductors contain only those edges 
corresponding to inductors and specified current sources. It ■ 
follows, therefore, that the fundamental cutset ec[uations 
introduce only zero entries in £522 ^23' ^22 

A 22 - it can be concluded that 0 ^^ = 0 and £52^ 

Likewise, since the cutsets defined by the edges of do not 
include edges corresponding to current variables i^'X|Jp 2 _/ 
follows that ^32 ~ 
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The state equations for the system are developed by 
substituting Eq.(A.6) into Eqs.(A,4) and (A, 5) with = 0^ 
^23 ~ ^0 9 - The results are 


f.?Pl 

dt 


~^1^^12 



°1 ^ 11%1 


“ ^1 ^10 ^pO 

and 

"Wp2"' 


u 


0 


U + D3 O33 


(A. 7) 


5= — 

r" — 

^2 ^^21^'^'^'^^^ 

PP3 - 

r~ 

D2 020 


°3 ®31 

L J 

1 

i 

^ 3^30 

L 


^ po 


(ATS) 


If the coEfficient matrix to the left of the equality 
sign in Eq.(A,8) is nonsingular, then Eq.(A,8) can be 
solved explicitly for 


^ p2 


^2 ^^21*''^'^'^^^ 

^2 ^20^'^^'^^^ 

1 ' 

r ~i 

i^pil 

^p3 


[ 

^^22 


"^pO 


p2 


and 




The result is 


(A. 9) 
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The above equations can be reduced to the follov/lng form 


diVpi 

dt 

= ^ h'^po it^po 

(A. 12) 

where 

p = [U - D;L ^12 °2 



, , -1 

Ql “ - °1 ^12 °2 >^21] "2 



-1 

□2 - [u - j 2)^2 ^2 ^ 21 ^ ^3 



if the coefficient matrix on the vector (d/dt)l|-)^j^ is non- 
singular. These coefficient matrices are indeed non-singular[l2]7 



An Example fo r I llus trat Ion s 


Consider the network shown in Fig. A.l. The system 
graph is shown in Fig. A, 2, The tree is selected maximally 
putting all the capacitors in the tree and all the inductors 
in the link',' 

■qjpo = 

"fpso “ Llo 


where i^ and i are the currents through the voltage sources 

V and V . 
o 


-Tpi 


T 

V V i i 1 ] 

3 6,7 8 10 





^6 ''7 ''8 '' l / 



1/Sc ° ' ° 

0 1/Cy,, 0 

O 0 1/C 

0 0 0 

0 0 0 

0 0 0 


0 

0 

0 

1/2L^ 

0 

0 


0 

0 

0 

0 

1/L^ 

0 



There are no capacitors in the cotree, ho inductors in the tr'^o 


Therefore 


Alp2 =° 



0 


(l)Vo 


(7) 

■i-Lt 

'^7 

C’ac 

a? 


CD 


i 




! '■ 

^ h- 

Ljvwv — 


O 


,-r?5?i?Qr\- 


— ^ Li 


><^i 


l__yv\-vv_i 


f 


'■^1 


■■ yc 

C3i 




CM.).4a ( 5 ) k ) V 

re 

"> K 

T A'' 

c 


- <^5 


0 


F\ 6} A ' 1 N £■ T {a)OR K wAiSi Ot. p 
(iQUiVALLMT C^^Ci^'iT U ^ IC 0 


i- or? 


1 fM 


Ti AJCt) 

ru K c N/ 


^TArt ^ 

>^yoL.TAc„ 



13 6 


~ QU' A T lc*ui 
iT C'''LCuL'P.T U." 


1 
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L^4 

^9 

’^ll 

■^33 ' ^4 

< 

VO 

^11- 



The circuit and cutset equations of the system graph are given 
by Hq.,(A.6) , '’For this exarrple, the nonzero matrices of 0 are 
as follows : 




13 8 


0 10 
1-10 



1 ° ±1 



— 

0 

-1 

1 

0 

0 

0 


0 

0 

1 

^31 

-1 

1 

0 

0 

0 


; j2! 

33 

0 

0 

0 


0 

0 

0 

0 ' 

0 

Jj 


-1 

0 

0 


Substituting the above defined matrices in Eq.(A.12) and 
simplifying ^the state equations are obtained v/hich are given 
by Eq. (A.13 ) . 
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APPENDIX B 

TIME DOMAIN SENSITIVITY CALCULATIONS [ b] 


The adjoint nc-tx^rork method will be used for sensitivity 

calculations in the time domain for linear dynamic netv'/orks. 

The problem to be investigated xin this section is the 

determination of VQ(t^)/c)Xj or OiQ(t^)/0Xj# where is 

some element parameter (R#L/C^gj^/ temperature, etc.)/ and 

V (t^) or i (t_) is the output voltage or currint at the 
o f of 

specific time t = t^T According to Tellegen's theorem, 

v^(t) 'i(r) = i^(t) v(r) = i^^(Dv(t) = v^(i:)i(t) = o (b.d 


where v(t)and i(t) are the element voltage and current vectors 
corresponding to the oi'iginal network (N), vC^) and '^('C) 
voltage and current vectors corresponding tos adjoint network 
The following equations can also be derived using 
Tollegen*s theorem. 


1"(T) A v(t) = 0 
v^(r-) Ai(t) = 0 
^■^(nAvCt) lv^(OAi(t) =0 


(B.2) 


where j!^\v{t) andAi(t) are the deviations of voltage and current 
of the perturbed network from the original network. 

Let i, (i +/51), and "i be the cerrents associated with 
the original network N, perturbed network Np, and adjoint 
Network^ respectively^ and similarly for the voltage vectors. 
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Subscript p Indicates port variables, and subscript b 
indicates variables for nonsource branches. With such 
partitioning of voltage and current vectors, Eq,<B,2)becomes 

i^(r)Av(t) - v^(r)Ai(t) = [rp(?)AVp(t) + i^(r) Av^(t)] 

**rVp(r) Aip(t) + 'v^(2:) Alj^(t)] 

= 0 

It follows that 

[-[ijr)fi5/p(t) - Vp(r ) Aip(t)] = [i^(r) AVj^(t) Aij^(t)] 

(B-.3) 

These equations are valid regardless of the values of t 
and ^ , In particular, Eq. (B.3) is valid if we let t and 
T. be related by ^ = t^-t, where t^ > 0 is the time instant 
at which v^ or i^ is to be investigated. With t and Z so 
related, further integration of Ect,(b, 3) on both sides 
from t = 0 to t = t^ leads to the following equation 

i ^[41 (T)Av (fe) + ^i_(t)] , dt, 

o ^ P P t^-t 

= / ^ [i^(i:)Av^(t) - V^(t) Aij^(t)]dt 

° C=t,~t 

(B.4) 

To evaluate the integrals in Eq.(B«4), first the 
transient analysis of N and N are performed to obtain the 


V 
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required voltages and currents. In the analysis, of N, the 

time variable is designated at t; in the analysis of’”^ 

the time variable is designated as XL • Responses of N are 

A. 

obtained for the Interval 0 t £ t^* the response of N 

; 

are obtained for 0 £ "C < Since C = tf-t/ as t increases 

from 0 tot: t^# the variable actually decreases from t^ to 

O',' For this reason# the adjoint netv;ork is sometimes said 
to be analysed in backward time# whereas the origin-^ 1 
network is analysed in forv/ard time. 


Let there be b non-source branches in N (and hence 

in The left side of Eq.(B,4) which is considered with 

the independent sources in N and N is evaluated as follows. 

It is considered any output voltage in N as the voltage 

across an independent current source# and any output current 

in N is considered to be the current through an independent 

voltage source'.' To get the left side of Eq,(B.4) as one 

term, which is precisely the desired change in output 

( ^v or Ai ) at t == t„. If the desired change to be Av , 
o o f o 

/\ 

the unit impuse function 6(^) is used and let 1^(D=*^(^) 
is applied across the output terminals and all other 
independent sources in ^ be zero valued. The L.H.S. of 
Eq.(B,4) becomes 


t ' f 

I ^ [-lp(?)Av_(t) + v^(r) Aip(t)] dt , = -I [i' (C)AV (t)] dt 

o ^ ^ ^ P O ° ° 

= / 6(t-t^)Av^{t) dt =Av^(t^). 
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If the desired change to be£^i^, the unit impulse 

function 6(r) is used and let'v^Cr) = 6(r) and 'all other 

A 

independent sources in N be zero valued. The L.h.S, 
of Eq, (B.4) becomes 


+ Vp(r) Aip(t)] dt^ =/ [v^(^)AiQ^t)] dt 


T=t£-t 




= 6(t-t^')Ai^^(t)dt =Ai.Q(t£). 


The R.H.S. of Eq',(B.4) in scalar notations can be written as 
t) t 

>?.1 i. , (B.5) 

° ? = tj-t 

Thus, the integrals in Eq.(B,5) are evaluated for each 
individual RLC branch and for each group of coupled branches. 
The results are then added to give the R.H.S. of Eq, (b., 5)V 
The evaluation of integrals for R and C is given below. 

For other types, the final results are shown in Table B.IV 


Resistance branch: 


v(t) = R i(t) and v((;:) = R ^(C) 


/ i(C)Av(t) -^CDAKt)] dt 
° r= t.-t 

tf ■ ^ ^ 

= / [x (C)[RAi(t) - i(t)AR] - R i(r)Ai(t)] dt 


C = t£-t 


[/ [i(C) i(t)] dt R 

® P" SS t “t 

^ f 


(B.6) 
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Capacitance branch 


l(t) = dq/dt, q(t) = Cv(t); and KZT) = d^dZI^'qCt) = c'vCO 

tjr 

A 


/ [irO^v(t) “ v(?) Zi»i(t)] dt 


Z. =t^-t 


= (T)^v(t) - v(2r) (iidq/dt)] dt 


C-t£-t 


= / pA . ‘■A 


[±(r) A.v(t)-v(^l) (d/dt )[CA.v(t) +v(t)AC j] dt 


C =t^-t 


(B.7) 


A 




The evaluation of of Eq.(B.7) is straight- forward. 

^ [v(r) v(t)Aic] dt 
° 2 =t .-t 

• 

= [- / [v^O v(t)] dt ]AC (B.8) 

° C =tf-t 

The evaluation of l3_+l2 Eq»(B.7) makes use of the method 

of integration by parts ^ which states 

b ,b . b 

X x(t) (dy/dt)dt = -/ (dx/dt)y (t)dt + [x(t)y(t) 

a a 

(B.9) 


we have 

I 1 +I 2 =/^ [i'(r)av(t)-vU) (d/dt[CAv(t>]] dt ^ 

° Z 

in evaluating we make use of Eq.(B-9) with x(t) = v(Z.) 

= v(t£“t), and y(t) = CAv(t).' Then, recalling that dt^-dc^ 


we have 
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Il+l2 


f 

= / [i(r)<av(t) “(dv/dr)c,!iv(t)] dt . 

o trt.-t 


~ [ (v(t')c^ v(t) ) • ^ ] 

o 

= / ^[l(C)/i!,v(t) ’ -'i(C)Av(t)] dt 

° ^ =t£-t 


c v(o)Av(t£) + C ''v(t^)£iV(o) 


A, 


= -C ^(o)^v(t^) + C V(t _) /ii.v(o) 


A 


(B.IO) 


Since nothing has been said about the initial condition 

A 

of a capacitor in the adjoint netvjork N/ we can exploit 

this freedom to achieve some simplification of the analysis 

A /V 

of N. We shall afgume that any capacitor in N always has 

a zero initial condition; i.e., 

v^(r) =0/ when ^ = 0 

Then Eq . ( B . 10 ) become s 

II +I 2 = c'v(t^)Zw(o) (B.ll) 

In many situations the initial voltage of the capacitor in 
N remains fixed when other element parameters vary • Then 
Av(o) = 0,and Eg.(B.7) finally becomes 


/ ['i(C’)‘^ v(t) - v(C)iAi(t)] dt 
° ^=t£-t 

= [- / ^ [v{x) v(t)] dt ' ]AC 
° ■^=t£-t 


(B.12) 
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The various steps Involved in calculating^Q(t^)/c)Xj 


orOi^Ct _)/c>x . are as follows: 


Step 1: 


Step 2 : 


Step 3 : 


Step 4 : 


Perform a transient analysis of the original network 
(N) for the time interval t = [O# t^] . Obtain 
i(t) or v(t) for resistive branches^ v(t) for 
capacitive branches^ and i(t) for inductive branches. 
In the adjoint network (N) , impose the following 
two conditions: 

1, All capacitor voltages and siinductor currents 

A. 

in N are zero at ^ = C* 

2', All independent sources iniN are set to zero 
except 'for = ~6(r) (impulse function)/, 

when (t^) is of interest, or'v^Cc:) = 6iz) 

whenAiQ^t^) is of interest. 

Perform a transient analysis of the adjoint network 

/V A 

(N) for the time interval [0,tjp]* Obtain iC'C) 
or'v(t) for resistive branches, for capacitive 

branches and'^('C) for inductive branches. 

Evaluate Eq.(B.4) to obtain the sensitivity compo- 
nents. With the excitations for ^ chosen as 
described in step 3, the L.H.S. of Eg.(B.4) should 
be exactly ^VQ(t£) onC^o(t^)'.' The R.H.S. is 
evaluated with /aid of Table B#l* 


step 5: Pindav^(t£)/axj or ai^(t£)/0Xj from the result of 
step 4, using the fact that df/c>Xj where x 

is the only parameter that varies, and^iXj 0* 
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Table B.'l Time Domain Sensitivity Components •' 


Element Description 
type 




/^[i{0^v(t) - v(r)Ai(t)] dt 

o 


^=tf~t 


R 


v_ = R i„ 
R R 


i^==dq/ dt 
q = Cv^ 


[/^{i^CD ij^(t)] dt ]ar 

° r=tf-t 

(r) V (t)] dt . >C 
o ^ ^ Z=t^-t 


G 


l=Gv^ 
G G 


o 


=t 


2:=tf 


Vt =d>/dt 
J-i 

)s= Li, 




yt^ 


V2=jAVi 
^1= 0 


.A 


[ / [i2(fl^) dt , ]A|A 

o cr=t£~t 


^2 ^^1 
''1=0 




Z=tf-t 


"m 


i_=g v- 
2 ^m 1 

vi=0 




T=tf-t 


m 


^ = 


r x, 
m 1 


V. 


= 0 


= nv 


1 

ni. 


A 


[ / [i2(C)iT_(t)] dt 


o 
t. , 


Z =t ~t 
r 


[ / [i2^2’) v^it) + v^(t) i2(t)]dt 


O 
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APPENDIX C 

EQUIVALENT CIRCUITS AND STATE EQUATIONS 
FOR DETAILED ANALYSIS OF DAMPING CIRCUIT 
LOSS CALCULATION";' 


The simplified equivalent circuits of the converter 
system during the sections 2 to 9 of Table 5;i are given 
in Figs', C.l to C.8. The state equations and the 
expression for current through valve 3 damping resistance 
for the section 2 are given by Eqs-(C.l) and (C.2) respec- 
tively, Similarly the equations for the section 3 are 
given by Eqs.(C.3) and^C.4) respectively. 








1^14 ' 

P^16 “ + Vg + 

PT'is “ '■'^18 ''7 ''s ''11 ''12 

P^20 “ ^”^20 '^7 '^8 '^11 

P ^22 = <^2 + ''7 ''s + '^11 + 

Pi4 ” 





(9 y'- 

^6> 


■7a 72. : Gd 


-fi^ G! S(^P^/)>^ fQouiScJ^vJ^ Cl^C^- ^he 

OGnOeif-ier Ccj^-^-e ^ ofa/c^^ % pe-ho^f ol-hA4 < cot <' "< 

C ISoJot.^ ^a<W '3 Ode Coior^tn/ €oi^- ) - 


'Ba Lt 


■ob Li 





S) ^ ® 

— -L (? ^o- 

- T^d -67 y-JT "■''*^1 

4 -^ ^ 4 . 




fxjCz . GcrnpU l^>')?_o p 010^0 . Ua,U- Qd CjuJ- <^ A.e c^(S>eitcyi 

ci.odPi-.. iho pc'^oS ^L^<s< KcoiCA'L-i-oi-f-/ 


Qfpe-i 


cJlt'i^^ ^l^e< :<Co£ <C A'l^-i-oi+A^ 

A , S L%.to <7^ LC^cUaJ :--'-^) . 



■^2 ^ Bc-mpliytd CxycjZJt- q />re Ccm 

cL()ou^ //?? />PT^oef < 0 ± < 2 ^/ 3 +-^ 

^ C3autdl6 OcJliM oxj> ■ 


Oa Lt 




151 



C 5 Sc^p ^Cpuu <S>oJevl' CaCuiJ- fh<f CouonhY Sy\ie.yy) 

Ciu^^Q /^>7 0C/ 2^/^ o(-i- -M < < A-4- 

pCCind I lQisUt<^ oX<? c/c^J) , 



A y C'Q SiYypli-p'<zJ Sq^iut o>af^xj- c'xUxxf- CanUi^ey S^j -^txM. 

cfiXtiO^ tfk p^ YJOcJ :T -f- oC < Cob ^ X ■ 

(^ G ^ ^ <r;uAx# 2. L^.A COti ) 
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- {l/3L^)[-v, - vg - + 2v^ - 


R 


'"'14 + "ll + ''12>/Pd 


(C.2) 


V7here 


V7 = [Q/(Q^-pr)]x ~ [p/(q^-pr)]y 
V ii = [r/(pr-q^)]x - [q/(pr~q^)]y 

and Y 

where P/Q/R^X/are defined as beloW; 


P = ( ~ + — ) 
^d 


2 2 

R= 

^d 


X - 1 + ^14 ~ ^12 ^ ^16 ~ ^8 ~ ^12 ^18 " ^8 ^ ^12 

' " ^ % ^G ^d 

. ^20 " ^8 “ ^12 . "^22 ” "^8 “ ^12 . ^ 


V •• 'V V - •- V '— V V mmm \r *mm i;r 

. 10 8 ^ 16 ^8 12 . 18 8 12 

4 R4 R^ R^ 


. ■ ^20 “ ^8 ~ ^^12 . "^22 “ ^8 “* ^12 1 ^ 

+ + 


The initial state variable vector for Eq.(C.l)is given 
by 


vg(t^) = -3vy2 
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v^g(ti) -3v^/2 






~3v^/2 




= “^d 


P^12 = 


pvi4 = ^-^14 + ■^11 + 

P^16 " ^“"^16 ^11 

pvi8 = i~^is + ^11 

1^20 = '-'^20 + '^11 ■" 

I"22 “ '■"'22 ■* "'ll ■*■ ''l2'''S''a 


(C.3) 
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where 

and 


pi 4 = (l/3L^)[-2v^,- - V,, + + vj 


a -11 '12 ■ "c ^ 'bJ 


pig = (l/3Lg)[-v^ 1 2v^^ - 2 v^2 - 


= <-14 - ''11 - ''l2>/''d 


(C.4) 


''11 = , 

t o t 

V — V v^v v-^v 

. .14 ^12 . "^16 U2 . ^18 12 

M = l£. H r H =; + r — 

5 R^ R j 


^20 “ "^12 , "^22 ~ ^^12 
+ ■ — ' + 


R, 




N = 


^ -f i 
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appendix d 

DggiyATION AND METHOD OF .SOLUTTQM OP T.TAPUNOV MATRIX 
EQUATION AND EVALUATION OF PERFORMANCE INDEX [s] 

Derivation of the Matrix Liapunov Ecruatlon 

Considering the unforced system 

X = A X (D.l) 

For stability of the above system there must exist a 
Liapunov function, or a v function, which is positive 
definite . 

Considering a general quadratic V function 
• T 

V(x) = X p X (D.2) 

where P is any arbitrary positive definite matrix- The 
time derivation of V(x) is given by 

^ • T T 

V(x) = X PX + XP X (D.3) 

and substitution of x from Eq,(D,l) gives 

- T T T 

v(x) = x A Px + x PAx 

= x"^ (A^P + P A)x (D.4) 

which can be written as 

V(x) = -x"^ Q X (D.5) 


whore Q is given by 
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A-* 


P + P A = ~Q 


(D.6) 


with any arbitrarily chosen positive definite P matrix# 

Q and hence V are# in general #indefinite , On the other 
hand arbitrarily selecting a positive semidefinite matrix Q 
and solving Eq.( i? 6 ), we get P# v;hich is positive semi- 
definite for an asymptotically sta- bl« system. 


So lution of the matrix Liapunov Ecfuation 

The matrix Liapunov equation is a linear equation 
in terms of the elements of the matrix P which has to be 
determined. Since P is symmetric matrix only m(m+l)/2 
elements of the P matrix are independent where m is the 
order of the matrix. 

The method of solution of Liapunov equation used 
here is based on rewriting the equation as 

[s] £ = 2 (D.7) 

whore, the matrix S is of order N = m(m+l)/2. If [s] and 
2 are known then r can be determined. 


The vector=x and 2 defined as 


£ = l^^ll' ^12' •**" ^in' ^22' 
2 = ^ 12 ' 


*/ 


rr 


rn ' ' nn^ 

(D.8) 






rn 


°ln' ^22' •**' 


(D.9) 
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the elements of vector r can be expressed in general^ 
form as 


= P L <K, for T, = 1/2,..., m 

UK 

K “*" 1,2,... ,m 

The matrix S is defined by 

NR = (2m ~ i + 2)^i - l)/2 + (j ~ i + 1) 
NC = 2m - L+2){L-1)/2 + (K - L + 1) 

then 


(a) 

for i = J 




S(NR, NC) 

= 2 Q 

KL 

for L = i, L i K 



= 2a 

IK 

for K = i, L < K 



= 0 

for U 1 

(b) 

for i < j 




S(NR, NC) 

= a. . + a 

. . for L = i/ K = j 



^Kj 

for L= i/ K j 



^Ki 

for Us: j, K > b , 


' 

= 0 

other .conditions* 

Evaluation of performance index 



The performance 

* 

index 

defined as 

can 

03 — 

= / x" Q X dt ( D *10) 

o 

be obtained by substituting value of Q from Eq.(- d - 


T 


OD m 

/ X (A P + P A)x dt 


or 
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00 

V? ‘ T T 

^ A P X dt 

o 

. ® T 

- / X P A X dt 
o 

From Eq, ( D .1) 



and thus Eq‘. ( D '.ll) can be written as 

'VL = X dt - / x'^ P X dt 


o o 


, m 


00 

GO „ 

00 

-[x 

P x] 

+ 

/ X P X dt 

- ; 



o 

o 

o 

T 


T 



X P 

X *" 

X 

P X 


o . 

o 

oo 

oo 



ForsE any asymptotically stable system 
'\= 


T 


X 


0 . 


(D 11) 


P X dt 


so we get 


( D 12) 
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APPENDIX E 

EQUIVALENT CIRCUIT FOR CALCULATION OF TURN-OFF OVERVOKPAGES 

AND VALVE DAMPER DESIGN 


The state equations for Pig.E.l are derived uing 
the algorithm given in Appendix A. These equations 
are written in the form 

px=Fx+Gu (E*l) 

where x is the state variable vector with 11 state 
varicibles given by 


[x] = [vg, V,, Vj, V^g. Vjj, I3, I5, l^g] 


(E.2) 

F is a 11 X 11 matrix with the following non-zero elements 


F(l,l) = [-(1 + 


^\( L. ^ ' i- 1. 1 .1/ 

‘^a *^3 ’’■■"''a 

^ ^ab ^ca 


F(l*2) [^SDC^^ca^^ ~ ”” R+R= 


G "'d “d 

C 

c . bo , 




F(l,3) = (1 + ^)(1 /Rg)/[c^(1 + ^ + 

° ^ ^ab ^ca 


ca 


^bc . 


F(l,4) = 

F(1,5) = (1+ 

ca ca 


)] 
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r'tiL-1. £:QufVALG/vT Circuit A-r Calculatidk; 
Tu6'ru OFF cAitR l/OLT ATnGS AVD VA'.va 

£ Si Cr-i KJ 


AAA,'- 
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F(l,6) = 

ab ca 


F(l,7) = (1/R,)/[C , (1 + 


^ + ^)] 
d--ab'- • ^ 




F(l,8) = [ll/(R+R^)' / [C^^( 1 + ^ + ^ )] 


- T. C 

ab ca 


S^c^ A r / *1 « =bc . Sx., 


F(l,9) = -Cl+ (1+ - 2 S + ^)] 

^ca ab ca 


FC1,10)= (C^/C^a>/f<^ab'l 


'ab 


ea 

F(l,ll) = [r./(R+R.)]/[C , (1 + ^ ^ )] 

^ ^ ^ab ca 


r(2.1) = [<Cj^/Cab>' fe ■" ■ i ■ 




d "d 

)] 


Sdc , 


'^ato ^ca 


P(2,2) = [-(1 +^)(|,+ !;> - i, - (VCr+r^,))]/ 


ab ^d ^d 

r 

be 


C C 

C^ca'l " ^ ^ C 


+ ^)] 

ab ca 


Sc . Sc, 


F(2,3) = 


% 


F( 2 , 4 ) = (1 +^)( 2 /R^) ^ ^ 


'ab 


f( 2,5) = fe 


f( 2,6) =[(1 + ^)(2/Ra)]''[Cea‘^ + ^ + 


'ab 


®ab =ca 


(E,3) 
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F(2,7) = (l/R<j)/[Coa(l + — ■^>1 

^ab ^ca 

F(2,8) = [-1/(R+Ra)]/[C (1+ ^ + •^)] 

ab ca 


f(2,9) 

F(2-;iO) 
f( 2 -^ 11 ) 




+ ^)] 


ca 


= -d . ^)/ [C^,(l ^ ^ ^,] 


'ab 


ab 


ca 


= [R,/(H-J]/[c .(1 


ca 


^ab *^ca 


F(3,l) = 1/C^Rq 

F(3,3> = -1/C^Rq 

F^4,2) = 4/C^Rq 

F(4,4) = -4/C^Rq 

F(5,1) =‘l/C^Ra 
F{5,5) = -l/C^R^ 

F(6,2) == 4/C^R^ 

F(6,6) = -4/C^R^ 

F(7,1) = 1/C^R<3 

F(7,2) = l/C^R^j 

F(7,7) =-l/c^Rd 
F(8A) = -1/C(R+R^) 
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F(8,2) = -1/c(R+R3^) 

F(8,8) = -1/C(R+R^) 

F(8,11) = R^/C(R+R^) 

F{9A) = 2/3L^ 

F(9,2) = l/3\ 

F(10,1) = l/3L^ 

F(10,2) = 2/3L^ 

F(ll,l) = -R^/.L^(R+R^) 

F(11,2) = -R /L (R+R^) 

a. oL a. 

F(11,8) 

f( 11A1)= -RR^/L^CR+Ra) 

G is a 11x2 matrix with the following non-zero elements'^' 
G(9,1) = “2/3L^ 

G(9,2) = -l/3L^ (E.4) 

G(10,l)= -1/3L^ 

G(10,2)= -2/3L^ 
u is a vector given by 

u = (Vj^ "^2 1 (E.5) 


where 



V 

Q, 



(E.6) 
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where 

V = sin (60 1 . 4 - 90 ) 

^ 3 ^ 

V, = ^ s±n(6^^t - '30) (E.7) 

a ' 3 L 

V = E_ sin( 'Ot - 150) 

C *3 L 


In Chapters 
'^01 ^"02 


3 and 4, the constant voltage sources 
are considered which are given by 


^01 " ^ 1 ^"''^^ " 

'^02 “ '' 2 ^ = “•*• O 


(E.8) 
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